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Abstract. We consider the construction and classification of some new mathematical ob- 
jects, called ergodic spatial stationary processes, on locally compact Abelian groups, which 
provide a natural and very general setting for studying difi^raction and the famous inverse 
problems associated with it. In particular we can construct complete families of solutions to 
the inverse problem from any given pure point measure that is chosen to be the diffraction. 
In this case these processes can be classified by the dual of the group of relators based on 
the set of Bragg peaks, and this gives a solution to the homometry problem for pure point 
diffraction. 

An ergodic spatial stationary process consists of a measure theoretical dynamical system 
and a mapping linking it with the ambient space in which diffracting density is supposed 
to exist. After introducing these processes we study their general properties and link pure 
point diffraction to almost periodicity. 

Given a pure point measure we show how to construct from it and a given set of phases 
a corresponding ergodic spatial stationary process. In fact we do this in two separate ways, 
each of which sheds its own light on the nature of the problem. The first construction can 
be seen as an elaboration of the Halmos-von Neumann theorem, lifted from the domain of 
dynamical systems to that of stationary processes. The second is a Gelfand construction 
obtained by defining a suitable Banach algebra out of the putative eigenfunctions of the 
desired dynamics. 



Outline 

This paper is concerned with mathematics of diffraction. More specificahy we are interested 
in the famous inverse problem for diffraction: given something that is putatively the diffraction 
of something, what are all the somethings that could have produced this diffraction. 

Diffraction has been a mainstay in crystallography for almost a hundred yearsj^ With 
the proliferation of extraordinary new materials with varying degrees of order and disorder, 
the importance of diffraction in revealing internal structure continues to be central. The 
precision, complexity, and variety of modern diffraction images is striking, see for instance the 
recent review article [53]. Nonetheless, in spite of many advances, the fundamental question 
of diffraction, the inverse problem of deducing physical structure from diffraction images, 
remains as challenging as ever. 

The discovery of aperiodic tilings and quasicrystals revived interest in the mathematics of 
diffraction, particularly since the types of diffraction images generated by such structures - 
essentially or exactly pure point diffraction together with symmetries not occurring in ordinary 
crystals - had not been foreseen by mathematicians, crystallographers, or materials research 
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scientists!^ Starting with [19], the resulting new techniques created in order to understand 
mathematical diffraction have been considerable. Among them we mention the introduction 
of ideas from the harmonic analysis of translation bounded measures [313], ergodic dynamical 
systems [3 |39l 001 HI] , and stochastic point processes jlTl HOl [1] . 

The diffraction of a structure is defined as the Fourier transform of its autocorrelation. The 
latter is a positive-definite measure in the ambient space of the structure - typically , but 
mathematically any locally compact Abelian group G - and the diffraction is then a centrally 
symmetric positive measure in the corresponding Fourier reciprocal space - either ~ 
or generally G as the case may be. 

In this paper we examine three basic questions that arise from this theory. The first 
question asks whether every positive centrally symmetric measure is actually the diffraction 
of something. The second asks about the classification of the all various 'somethings' that 
have this given diffraction. The third asks what kinds of structures the 'somethings' can be. 
In the case of pure point diffraction (when the diffraction measure is a pure point measure) 
we can give a satisfactory answer to these problems. 

We will be in the setting of locally compact Abelian groups. The difficulties of the questions, 
at least as we examine them here, have little to do with the generality of the setting. They 
are just as hard for M^. Indeed the second question is a very old one in crystallography, called 
the homometry problem and is part of the fundamental problem of diffraction theory, namely 
how to unravel the nature of a structure that has created a given diffraction pattern. 

Our first question immediately raises the third question. Exactly what do we mean by 
a structure that can diffract? Typically diffractive structures are conceived of as discrete 
measures, for instance point measures describing the positions of atoms, vertices in tilings, 
etc., weighted by appropriate scattering strengths, or as continuous measures describing the 
(scattering) distribution of the material structure in space. More recently, starting with the 
work of Gouere |16| [T7] there has been a shift to discussing the diffraction of certain point 
processes, which in effect are random point measures, and various distortions of these (see 
e.g. [U [ini [29] ) • In this paper we are led to introduce new structures which we call a spatial 
stationary processes. We shall show that these always lead to autocorrelation and diffraction 
measures and contain the typical theories of diffraction discussed abov^ In terms of these we 
get an attractive answer our three questions, in the case of pure point diffraction. However, 
the answer to the third question is interesting, because it does not describe the underlying 
structure explicitly, but rather in terms of the ways that it appears to compactly supported 
continuous functions on the ambient space. In this respect it resembles the description of 
measures given by the Riesz representation theorem or the theory of distributions with their 
spaces of test functions. In our case, the convergence conditions of measures do not seem to 
be generally valid. In this respect there is much to be learned about the full scope of which 
structures can diffract. ^ 

Let G be a locally compact Abelian group and G its dual group. We treat these groups 
in additive notation. We begin with an abstract notion of a stationary process for G and 
then shift our attention to the more concrete notion of a spatial stationary process for G (§2^ 



and ^3.3). As for terminology let us note that outside of the discussion of usual stochastic 



^Indeed, the discovery of quasicrystals by Shechtman via diffraction experiments met with substantial dis- 
behef for a period of about two years. Ifis determination to communicate it to a skeptical scientific community 
is particularly stressed by the committee awarding him the Nobel prize in Chemistry in 2011 for this discovery. 

n 

There is a technical requirement about the existence of second moments. 
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processes, the word 'stochastic' is not used in our discussion of processes and stationary 
processes. 

A spatial stationary process basicaUy consists of a G-invariant mapping 

N:Cc{G)^L\X,fi) 

where {X, ji) is a probabiUty space on which there is a measure preserving action of G. 

The idea here derives directly from the theory of stochastic point processes, and it is perhaps 
useful to discuss this briefly. We follow |22j here. A stochastic point process, say in W^, is 
a random variable M on a probability space P) with values in the space Mp of point 

measures ^ on for which ^{B) e Z>o for all bounded subsets B of W^. These are point 
measures and can be written in the form ^ = X^j^i where K can be a positive integer or 
infinity and the points Xj are (not necessarily distinct) points of W^. Usually the space M.p is 
far larger than needed. For one thing it allows points in to appear with multiplicities and 
for another it does not require any minimal separation between distinct points (a hard core 
condition) that would usually be imposed in crystallography or in tiling theories. We assume 
then that we are really only interested in some subset X C Mp of permissible outcomes for the 
point process, and that almost surely the point process produces measures in X. Typically we 
wish X to be invariant under translation since the position of our point sets is not particularly 
relevant, and we assume that the law /i of the process (the probability measure /x induced on 
X by the random variable) is invariant under translation (i.e. the process is stationary). 

Now the point is this. If / G Cc(M'^) then we obtain a random variable on (17, J^, P) by 
^fiO ■■= Ej=i fixj), where C = M{u) G X is the measure indicated above. In effect N 
can be viewed as a mapping from Cc(M'^) into mappings on X. Under mild assumptions N : 
Cc(M'^) — >• L^(M'^,/i). As is often the case with random variables, we ignore the underlying 
probability space (il, P) since we have everything we need from the probability space {X, fi) 
and the mapping N. Then N itself can be referred to as the stochastic point process. 

In this paper M"' is replaced by any locally compact Abelian group G. Also X is not 
assumed to have any particular form (e.g. a set of point sets, or a set of measures). We do 
not wish to assume any particular mathematical structure on the outcomes of the processes 
we are studying. Our point of view is that we will start with a measure in G which we would 
like to show is the diffraction of something. But we do not know in advance what sort of 
mathematical objects could produce this diffraction (point sets, continuous distributions of 
density, Schwartz distributions, etc.). We derive information about X only from the behaviour 
of the functions N(F), F S Cc(G), upon it. This is, physically, rather appropriate. In a 
physical situation one determines the structure of a physical object by measurements upon it 
and the results of measurements are ultimately all that one can know about it. 

We are particularly interested in pure point diffraction. However, the concept of stationary 
processes seems interesting in its own right and we develop the theory of these a bit, indepen- 
dent of the question of pure point diffraction. It is even useful initially to ignore the fact that 
the resulting process is spatial, that is, it is supposed to result from some structure in G. So 
at the beginning we often only require that we have a G-invariant mapping N : Cc(G) — > % 
where % is a. Hilbert space with conjugation - we call these stationary processes. Assuming 
existence of second moments, we can then, for each such process, find an 'autocorrelation 
measure' on G and a 'diffraction measure' on G. The latter is always a centrally symmetric 
positive measure, and through it we can define continuous and pure point diffraction. We 
show how pure point and continuous diffraction are tied to notions of almost periodicity of 
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A^. Also we show how to decompose a sum of spatial processes into sub-processes which are 
respectively pure point and continuous. 

Coming to the pure point part of the paper, ^ we assume given a centrally symmetric 
positive pure point measure oj on G, which can be expressed as a Fourier transform. We wish 
to create a probability space {X, fj,) on which there is an ergodic action of G and a continuous 
stationary process N : Cc(G) — > L'^{X,fi) whose diffraction is u. 

Not surprisingly phase factors play a crucial role in this since it is phase information that is 
lost in the process of diffraction. There is a natural combinatorial type of group Z generated 
by the positions of the Bragg peaks. This group can be seen as an abeleanized homotopy 
group of an associated Cayley graph. It is the dual group of Z, or the dual group of a natural 
factor Z of this group, that classifies all the spatial processes on G with diffraction uj, the 
classification being up to isomorphism or up to translational isomorphism respectively. 

We call these phase factors phase forms. Background on phase forms is given in Section 
[7[ In ^10 we show how to use u and a phase form a* to construct a spatial process whose 



diffraction is lo. Uniqueness is already shown in ^ How all these ingredients give a solution 
to the homometry problem is discussed in 

Our solution to the pure point problem is in some sense an elaboration of the famous 
Halmos-von Neumann result about realizing ergodic pure point dynamical systems in terms 
of the character theory of compact Abelian groups. There the objective was to classify pure 
pointedness at the level of the spectrum. In our case we wish to classify the diffraction, which 
amounts to not only knowing the spectrum but also the intensity of diffraction at each point 
of the spectrum. Our approach to the construction of a spatial process is to realize X as 
the dual to the discrete group generated by the set of Bragg peaks, but at the same time to 
include all the phase information which derives from the Bragg intensities. In this context 
we mention the recent work of Robinson [45j on how to realize systems with a given group of 
eigenvalues via so called cut-and-project schemes. 



In { 12 we take a closer look at the situation that the underlying group is compact. In this 



case all processes can be realized as measure processes under some weak assumption. We then 



specialize even further and in ^ 13.3 discuss some results of Griinbaum and Moore on rational 
diffraction from one dimensional periodic structures. These pertain to the simple case that 
G = Z/NZ, but even so, the results are interesting and not easy to obtain. 



Finally in { 14 we sketch out a second approach to the construction of a spatial process, 



this time based on Gel'fand theory. 
Summarizing: 

• Each pure point ergodic spatial process gives rise to a pure point diffraction measure 
LO and a phase form a* . 

• Spatial processes with the same pure point diffraction measure and the same phase 
form are naturally isomorphic up to translation. 

• The set of all possible phase forms associated with a given pure point measure w form 
an Abelian group 2{uj), and for each choice of a* € 2{uj), there is a ergodic spatial 
process with diffraction oj and phase form a*. 

The paper has two main features, after the introduction of stationary and spatial processes: 
one is a development of a general theory of these types of processes and the other is a detailed 
study of pure point diffraction and examples of how it looks in special cases. Readers interested 
primarily in the pure point theory can skip sections 3.3, 3.4, 4, 5 on first reading if they so 
wish. 
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1. Notation 

In this section we gather some notation used throughout. 

N := {1,2,...,}. Let G be a locally compact abelian group. We let Cc(G) denote the 
space of compactly supported complex valued continuous functions on G. For compact K C 
G, Cft'(G) is the subspace of elements of Cc(G) whose support is in K. For F £ Cc 



F G Cc(G) is defined by F{x) = F(—x). The translation action of G on itself is denoted by 
r: r(t)(x) = t + x for all t,x G G. This action determines, in the usual way, an action on 
Cc(G). By /(G we denote the (fixed) Haar measure on G. 

The convolution of two functions F,G £ Cc(G) is defined to be the function F*G £ Cc(G) 
given by 

F*G{t)= [ F{t-s)Gis)dlG{s). 
Jg 

The dual group of G is denoted by G and the Fouriert ransform of an F G Cc(G) is denoted 
by F i.e. 

Hi) = [ hJ)Fit)dlG{t). 
Jg 

A measure on G is called transformable if there exists a measure cr on G with 

iFPdcr = / F*Fdiy 



for all F £ Cc(G). In this case a is uniquely determined and called the Fourier transform of 
V. In such a situation we call a backward transformable. Unlike all other notation introduced 
in this section the concept of backward transformability is not standard. However, it is clearly 
very useful in our context as it is exactly the property shared by the diffraction measures we 
investigate. 

A sequence {A^} of compact subsets of G is called a van Hove sequence if for every compact 

(1) hm 7 = 0. 

Here, for compact A, K, the "i^-boundary" A of A is defined as 

d^A := {{A + K) \ A°) U (Gr\A -K)r^A. 

The existence of van Hove sequences for all a-compact locally compact Abelian groups is 
shown in [361 p. 249], see also Section 3.3 and Theorem (3.L) of ^Sli Appendix]. 

2. Stationary processes introduced 

In this section we introduce our concept of stationary process on the locally compact 
Abelian group G. To do so we first discuss some background on unitary representations. The 
concept of stationary process is somewhat more general than needed in the remainder of the 
paper. The reason for this is twofold: On the one hand this general treatment does not lead 
to more complicated proofs but rather makes proofs more transparent. On the other hand 
for a the study of mixed spectra this general concept may prove to be especially useful. 

Let G be a locally compact abelian group and let T-Lhe a. Hilbert space with inner product 
(•|-). We assume linearity in the first variable and conjugate linearity in the second. A 
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strongly continuous unitary representation of the group G on is a map T from G into the 
bounded operators on V. such that r i— {Trf,Trf) is constant on G for any / G "H, Tq is the 
identity on Ti, Tt+s = TtTg holds for all s,i G G and t ^ Ttf is continuous for any f £ 71. 
Then, obviously each Tg is unitary (i.e. bijective and norm preserving). An f € T-L is then 
called an eigenfunction of T with eigenvalue ^ G G if Ttf = {(,,t)f for every t G G. The 
closed subspace of Ti spanned by all eigenfunctions is denoted by Tip. The representation T 
has pure point spectrum if Ti = Tip, or equivalently if Ti has an orthonormal basis consisting 
of eigenfunctions |^ 

By Stone's theorem, compare [331 Sec. 36D], there exists a projection- valued measure 
Et ■ Borel sets of G — )■ projections on T-l 

with 

(2) {f\Ttf) = /(e,t)dp/e) 

Jg 

for all t G G and f £ Ti, where pjr is the measure on G defined by Pj{B) := (/ | ExiB)/). 
Then is called the spectral measure of /. It is the unique measure on G satisfying ([2]). 
A mapping (•) of a Hilbert space {Ti, {■ \ •)) into itself is said to be a conjugation of T-L if 

(i) (•) is conjugate- linear and ^ = ^; 
(fi) for ah en, = 0- 

Definition 2.1. Let G be a locally compact abelian group. A stationary process or T-L- 
stationary process on G is a triple M = {N,1-L,T) consisting of a Hilbert space % with 
conjugation, a strongly continuous unitary representation T of G on and a linear G-map 
iV : Cc{G) — > n such that for ah F G Cc(G) N(F) = N{F). A stationary process is cahed 
ergodic if the eigenspace of T for the eigenvalue is one-dimensional. 

The set of all ?^-stationary processes on G forms a real vector space in the obvious way 
by taking linear combinations of mappings. There is also a canonical notion of isomorphism 
between stochastic processes: The processes {Ni,'Hi,Ti) and {N2,'H2,T2) are called isomor- 
phic if there exists an invertible unitary map U : Hi — > T-L2 which intertwines Ti and T2 and 
satisfies N2 = UNi. 

We will be interested in the spectral theory of stationary processes. An important notion 
is given next. 

Definition 2.2. A stationary process J\f = {N,T-L,T) is said to have pure point spectrum if 
the representation T of G on % has pure point spectrum. 



Definition 2.3. A stationary process N = {N,1-L,T) is called a spatial stationary process if 
there is a probability space (AT, p) and an action T : G x A — > X, {t,x) t ■ x such that p 
is invariant under the G-action on X through T,'H = L'^{X,p) with the usual inner product 

■^Often the eigenfunction condition is written as Ttf = {^,t)f for eigenvalue £, £ G. For our purposes things 
work out more smoothly by formulating it with the conjugate value. 

^The origins of the definition lie in the theory of stochastic point processes, hence the name. We will often 
simply write process instead of stationary process. We do not use the word stochastic, although our definition 
allows for processes that are stochastic in nature. 
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if \ 9) = fx action of T is extended to an action on V. by Ttf{-) = /((— 

for all i G G, / G L'^{X,fi). The conjugation here is the natural one: complex conjugation of 
functions on X. We write M = {N, X, ^,T). The stationary spatial process is called full if 
the algebra generated by functions of the form tp o N{F), G Cc(C), F G Cc(G), is dense in 

L\x,^l). 

Remark 2.4. In some sense the assumption of fullness is only a convention. More precisely, 
given any spatial stationary process we may create a full process with essentially the same 
properties by a variant of Gelfand construction. Details are given below in Theorem |5.4[ 

For the sake of economy, we have used the same notation, T, for the action of G on X and 
the corresponding action of G on := /x). We often omit the word 'stationary' in the 

sequel, but we will always understand it to be in force. 

The definition of stationary processes implicitly defines them as real processes: N{F) is 
real for all real- valued functions F on G in the sense that N{F) = N{F). In the case of 
spatial processes, all real- valued functions on G are mapped by N to real- valued functions 
on X. One could relax the conditions to allow complex- valued processes, but we shall not do 
that here. 

Definition 2.5. Two stationary spatial processes M = {N, X, /x, T) and M' = {N' , X' , fi' , T') 
are spatially isomorphic if there is an invertible unitary mapping M : L'^{X, fj.) — > L^(X', /i') 
with 

• MoN = N'; 

• MoTt = Tl o M for all t G G; 

• Mifg) = M{f)M{g) for aU /,5 G and M-^/V) = M-\f')M-\g') for 
all f',g'GL^iX',^i'). 

The map M is then called spatial isomorphism between the processes. 

The isomorphism class of = (N, X, /x, T) is denoted by [M] = [{N, X, /i, T)]. 

Remark 2.6. (a) The first two points of the definition just say that the processes are isomor- 
phic. The special requirement for spatial isomorphy is the third point. As shown in the next 
proposition, the spatial isomorphism are exactly those isomorphisms which are compatible 
with the L°° module structure of in the sense that M{fg) = M{f)M{g) for aU f £ L°° 
and g £ L^. 

(b) Our definition of spatial isomorphism parallels the notion of spectral isomorphism for 
measure preserving group actions on probability spaces. Our conditions on M imply that 
spatial isomorphism implies conjugacy in the sense that there is a measure isomorphism 
between the Borel sets of X and those of X' [52], Theorem 2.4. If in addition X,X' are 
complete separable metric spaces then this can be improved to being an isomorphism between 
X and X' , that is a bijective intertwining mapping between two subsets of full measure in X 
and X' respectively [52J, Theorem 2.6. 

Proposition 2.7. Let Af = {N,X,fx,T) and M' = {N' , X' , fi' ,T') be spatially isomorphic 
processes with spatial isomorphy M. Then, M maps L°°{X,fi) into L°° {X' , and M{fg) = 
M{f )M{g) holds for all f G L^{X, fj.) and g G L'^{X, fi). The corresponding statements hold 
for M"^ as well. 
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Proof. It suffices to consider the statements on M. By assumption we have that for f,g£ L°° 
the product M{f)M{g) belongs to and in fact equals M{fg) (and similarly for M~^). By 
a simple limit argument, we then see that for any / £ L°° the equation M{fg) = M{f)M{g) 
must hold for any g £ L^. This shows that the operator of multiplication by M{f) is defined 
on the whole of (as M is onto). Furthermore, the graph of every multiplication operator 
can easily be seen to be closed. Thus, the closed graph theorem gives that the operator of 
multiplication by M(/) is a bounded operator on L^. This in turn yields that M{f) belongs 
to L°°. □ 

3. Diffraction theory for stochastic processes 

In this section we show how each stationary process with a second moment comes with a 
positive definite measure 7 on G called the autocorrelation, a positive measure u onG called 
the diffraction measure, and the diffraction-to-dynamics map 9. We then discuss how the 
autocorrelation comes about as a limit and characterize existence of a second moment. 



3.1. Autocorrelation and diffraction. Let a stationary process M = {N,T-L,T) be given. 
We say that M has a second moment (or is with second moment) if there is a measure 
:= iJ,^^^ on Cf-{G x G) (necessarily unique if it exists) satisfying 

(3) fi^^\F0G) = {N{F)\N{G)) 

for all real valued F,G £ Cc(G). Much of this paper is based on the assumption that the 
second moment measure exists. Since N is real, so is /U^^-*, and hence it is real valued]^ 
Furthermore evidently translation bounded since is G-invariant. 

Proposition 3.1. Let N = (NjT-LjT) be a process. Then N has a second moment if and 
only if there is a measure 7 = 7^^ on G satisfying 

(4) -f{F*G) = {N{F)\ N{G)) 

for all F,G G Cc(G). In the case that such a measure exists, it is unique and furthermore 
it is real (i.e. assigns real values to real valued functions), positive definite, and translation 
bounded. 

Proof. It suffices to restrict attention to real-valued F and G in Cc(G) as N is invariant under 
a conjugation and hence both F * G and {N{F) \ N{G)) are linear in F and antilinear in G. 
Let ip := /U*^^-* above. The assumption of stationarity of N shows that ip is invariant under 
simultaneous translation of both of its two variables. Define the mapping 

(5) G X G — > GxG 

{x,y) ^ {x -y,y) =: {u,v) 

For Fi,F2 G Gf{G) we define H G Cc(G x G) by H{u,v) = H{x -y,y) = Fi{x)F2{y). This 
change of variables defines a new measure ^* via 



Fi{x)F2{y)di,{x,y)= / / H{u,v) {u,v) 



^The mth moment of A/" will be defined by (Fi ® • ■ • ® F™^) = iV(Fi) ■ ■■N{F^)d^ which is difi'erent 
from what we get from the inner product (• ] ■) when m — 2 unless we restrict to real- valued functions. This 
is why we restrict to real-valued functions in this definition. 
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Since {Ttx,Tty) -f-T- {u,Ttv), the translation invariance of translates to translation invari- 
ance of the second variable of ip* . Thus for measurable sets A, B G G we have ip*{A x B) = 
'ip*{A X (t + B)) for all t G G. For A fixed this is leads to a translational invariant measure 
on G which is hence a multiple c{A) of Haar measure /(g: 

r{A X B) = c{A)Ig{B) . 
The mapping A i— )• c{A) is a measure on G, and this measure is the reduction ■tp^^'^ of Tp: 

r{AxB)=i;''''^{A) Ig(B). 
Now for H{u,v) = Fi{u + v) F2{v) , 

H{u,v)di;*iu,v) = j j Fi{u + v)F2{v)dip''''^{u)dlG{v) 

Fi{u + v)F2{v)AIg{v)^ dij'^^'^{u) 

(Fi*F2)(n)dV'«d(n), 
so finally 

and -0^^*^ is the desired measure, which we now denote by 7. 

It is positive definite since 7(F * F) = iP\F ® F) = \\{N{F)\\'^ > 0. Ah positive definite 
measures are translation bounded [8]. Since ip := fi^"^^ is a real measure, one sees that 7 is 
also a real measure. The denseness of the set of functions F * G in Cc{G) shows that 7 is 
unique. 

We can reverse all the steps of this proof, and from the existence of the real positive definite 
measure 7 satisfying Q derive the corresponding second moment measure 

Existence of a second moment immediately implies some continuity properties of A^. 

Corollary 3.2. Let {N,T-L,T) he a stationary process with second moment. Let 1 < p,q < 00 

be given with 1/p + 1/q = 1. Then, for any compact K C G, there exists a Ck > with 

\{N{F) I N{G))\ < 

for all F,G £ Cc{G) with support contained in K. In particular the map N : Cc{G) — > 
L'^{X,m) is continuous (with respect to the L'^-norm on Cc{G)). Also N : Ck{G) — > 
L'^{X,m) is continuous with respect to the sup norm on Ck{G). 

Proof For F,G £ Cc{G) with support contained in K the support of F * G is contained in 
K-K and ||F*G||oo < ||i^lUp(G)l|G||L,(G). As \{N{F), N{G))\ = |7(F* G)|, this easily gives 
the first statement. Now, the second statement follows by taking p = q = 2. Finally, for 
F G Ck{G), ||F||i2(G) < ||F||oo(/gW)'/'. □ 

The continuity property of the preceding corollary allows one to extend the map N. This 
is discussed next. 



Corollary 3.3. Let M = {N,T-L,T) be a stationary process with second moment. The map N 
can be uniquely extended to the vector space of all measurable functions f on G all of which 
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vanish outside some compact set and are square integrable with respect to the Haar measure 
on G. This extension (again denoted by N) is G-equivariant and satisfies 

I F*Fd-f = {N{F)\N{F)) , 
Jg 

meaning in particular that the integral exists and is finite. 

Proof. It suffices to show that N can be uniquely extended to L'^{U, Ifj) for any open U in G 
with compact cfosure. Let such an U be given. 

Let F G L?'{U,ljj) be given. As, U is open, we can then find a sequence {Fn} in Cc{U) 
converging to F in the sense of ([/,/[/). By the previous corollary, we infer then that 
{N{Fn)} must be a Cauchy-sequence, whose limit does not depend on the choice of the 
approximating sequence {Fn\. This shows that N can be extended to L?'{U,lu). 

By construction of N{F) and the definition of 7 we have furthermore 



{N{F) I N{F)) = lim / Fn * Fndj. 



Moreover, a direct application of Cauchy-Schwartz inequality shows that Fn * Fn converges to 
F * F with respect to the supremum norm and has support contained in U — U. This easily 
gives that 

lim Fn* Fndj = F* Fd-y. 
Jg Jg 

This finishes the proof. □ 

Remark 3.4. Note that any bounded measurable function with compact support belongs to 
Lp'iU, dl^) for any open U containing the support. Thus, N can in particular be extended to 
bounded measurable functions with compact support. 

The above discussion shows that any stationary process with a second moment gives rise 
to a real positive definite measure 7. As 7 is positive definite it is transformable i.e. its 
Fourier transform oj exists [HJ [50], and since 7 is a real measure, uo is centrally symmetric 
{oo{A) = ijj[—A) for all measurable sets A). The measures 7 and uj lie at the heart of our 
investigation. 

Definition 3.5. Let {N,T-L,T) be a stationary process with second moment. Then, 7 := ■~f]\f 
is called the autocorrelation of the stationary process. Its Fourier transform uj := wat := tjv 
is called the diffraction or diffraction measure of the stationary process. 

Definition 3.6. A stationary process {N,T-L,T) with second moment is said to have pure 
point diffraction (resp. continuous diffraction) if the diffraction measure is a pure point 
measure (resp. continuous measure). 

The definition of oo yields that it is a centrally symmetric positive measure and is the unique 
measure satisfying 

(6) / |F|2da;= / F*Fd7, 

JG Jg 

or equivalently, by linearizing this, 

' FGduj = [ F*Gd-f 
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for all F,G Cc(G). Considering the Hilbert space L^(G,a;) with the corresponding inner 
product being written (• | •)g and using the definition of 7, we obtain for any F,G €z Cc(G), 

(7) (F|G)g= [ FGduj= [ F*Gd-f = {N{F)\N{G)). 

Jg Jg 

This shows that the map F 1— t- N(F) is an isometry. We define an action T of G on L^(G, w) 
by defining Ttf for all / G L^(G,w) through 

ftfik) = V^fik) 

for all t G G and all /c G G. 

Remark 3.7. We note that if is a ?^-stationary process on G with second moment then 
for all c G M, cN is another such process. The corresponding autocorrelation and diffraction 
measures are then scaled by |cp. 

3.2. The diffraction to dynamics map. Any process with a second moment comes with 
an isometric G-map 6 from L^(G, io) to Ti. This allows one to transfer certain questions from 
n to L^{G,u)). Also, it means that for certain eigenvalues there are canonical eigenfunctions 
available. These topics are studied next. 

Proposition 3.8. Let N = {N,T-L,T) he a stationary process with second moment. Then 
{F : F G Cc(G)} is dense in L'^{G,(jj) and there exists a unique isometric embedding 

(8) e = 6n : L^{G,uj) 

with F I— )• N{F) for each F G Cc(G). This mapping is a G-map and furthermore, for all 
D G L2(G,a;), 

eiD) = e(D). 

Proof We begin with the denseness statement. Let L be the closure of {F : F G Cc(G)} 
in L'^{G,uj). We wish to show that L = L^(G,a;). As w is a translation bounded measure, 
the space Cc(G) is dense in L^(G,a;). Thus it suffices to show that Cc(G) belongs to L. 
As the convolutions of elements from Cc(G) belong to Cc(G) as well, we infer that for all 
Gi, . . . , Gm £ C'c(G) 

Fi ■ ■ ■ Fn ■ Gi ■ ■ ■ Gm £ L 

for all Fi G Cc(G). Now Fj belongs to the set Co(G) of continuous functions on G vanishing 
at infinity. An application of the Stone- Weierstrass theorem then shows that 

fGi ■ ■ ■ Gm £ L 

for each / G Co(G). Fixing / G Cc(G) we can then use another application of the Stone- 
Weierstrass theorem, to infer that 

fg^L 

for all / G Cc(G) and g G Co(G). This yields Cc(G) C L as we intended to show. 

The existence of 9 now follows by defining it initially hy 9 : F ^{^) for all F G Cc(G), 
using the fact that it is an isometric mapping by ([T]), and then extending it to the L'^-closure 
L which is L'^{G,lo). 
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As for the second statement, first we note that the G-action on L^(G,cj) is designed to 
make the mapping F F into a G-map, while F i— N{F) is already a G-map, and second, 

for all F G Cc(G), F = F, from which 



e \^Fj = e{F) = N{F) = N{F) = 9{F) . 

We finish using the denseness of the first part. □ 

Definition 3.9. The map 9 = On associated to a stationary process M = {N,T-L,T) with 
second moment is called the diffraction to dynamics map. 

The relevance of 9 for spectral theoretic considerations comes from the following conse- 
quence of the previous proposition. 

Theorem 3.10. Let N = {N,T-L,T) be a stationary process with second moment and 9 the 
associated diffraction to dynamics map. Then, for any H G L'^(G,lo) the spectral measure 



j.^^ is given by 



h'eiH) 

Proof. As ^ is a G-map and an isomorphism a short calculation gives 

{9{H) I Tt9iH)) = {0{H) I 9{{tr)H)) = fit,0\H\^doj. 

Jg 

Thus, the (inverse) Fourier transform of is given by t i— t- {9{H) \ TtO{H)). By the 

characterization of the spectral measure in ^ the theorem follows. □ 

From the previous theorem it is clear that if % has pure point spectrum then oj is discrete 
and the diffraction is pure point. However, it is not clear that pure point diffraction implies 
pure point spectrum but as we shall see, it is so: each implies the other in the case of spatial 
processes. 

The diffraction to dynamics map connects pure point diffraction to eigenvectors of the G 
action on %. Later, we shall see that in the case of spatial processes, this allows us to make 
a complete correspondence between pure point diffraction and pure point dynamics. 

For a pure point diffractive stationary processes with diffraction measure oj we define the 
set of its atoms by 

S = S{uj) := {A: G G : uj{{k]) / 0} = {A: G G : uj{{k]) > 0}. 

The set S is often called the Bragg spectrum of the process, and its elements are sometimes 
called Bragg peak^ 

In the sequel we will often omit the brackets when dealing with one element sets of the 
form {A;}. In particular, we will set U}{k) := a;({A;}) for A; G G. 

Let Ifc be the characteristic function of the set {A;} C G i.e. lfc(A;') is 1 or according as 
k' equals k or not. It is easy to see that these functions are the only possible eigenfunctions 
for our action of G on L^(G,a;). Define /fe. A; G 5, by = ^(Ifc). Then each is an 
eigenfunction in % for the eigenvalue k, in the sense that Ttfk = {k,t)fk. 



'''Sometimes, the term Bragg peak is also taken to mean both the position k of the atom and its intensity 



c.({fe}). 
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Lemma 3.11. Let N = {N,T-L,T) be a stationary process with second moment. Then, for all 
k e S, -k e S, and f^k = fk and \ \f-k\ \ = \ \fk\ \ = uj{ky/'^. 

Proof. Since u is centrally sy mme tric, u}{{—k}) = uj{{k}), which gives the first statement. 



For the second, using Prop. 3.8 we have = 6{l^k) = ^'(Ifc) = ^(Ifc) = fk- Then 



Wf.kW = if-k I f-k) = {fk I fk) = WfkW = {lkAk)=^{k). □ 

We can use the preceding considerations to compute the map N in the case of pure point 
diffraction. 

Proposition 3.12. Let {N^T-L^T) he a stationary process with pure point diffraction and 
associated Bragg peaks S. Let he the associated diffraction to dynamics map and fk = ^(Ifc)? 
k £ S. Then 

NiF) = Y,mfk 

kes 

for all F £Cc 



Proof. By definition of 9 we have N{F) = 0(F). Obviously, F = ^^gs -^(^)-'-fc and the claim 
follows. □ 

3.3. Spatial stationary processes: Two-point correlation. In this section we show how 
the autocorrelation can be given a meaning that agrees with 'classical' two-point correlation 
associated to stationary point processes. 

Assume that we are given an ergodic spatial stationary process J\f = (N, X, fj,,T). As dis- 



cussed in Corollary 3.3 we can assume that N is defined on all measurable bounded functions 
with compact support. We also assume that G has a countable basis of topology and is hence 
metrizable. Let be a van Hove sequence for G (as discussed in Section [T]). We 

assume that this is fixed once and for all. 

Definition 3.13. For G X we define the two-point correlation or autocorrelation of ^ at 
any F G Cc(G) as 

i™nrTm f ^1?^ TJJ-^ I iN{lB)N{F)){T^,Odk{x)] 
whenever the limit exists. 

This needs some comments: The limits in the definition are taken in the order indicated: 
first the inner and then the outer. B is an open (or measurable) neighbourhood of in 
G. The statement i? — )• means that we take a nested descending sequence {-Bm} of such 
neighbourhoods, all within some fixed compact set and that l<Q{Brn) 0. We are using 
the notation dx to stand for the longer d?(G(x). The definition requires that we give meaning 
to N{\b) as a measurable function on G. This us es th e extension of N from Cc(G) to 
L^-functions with compact support given in Corollary 



3.3 



The intuition behind the definition is as follows. The two-point correlation at ^ E X for 
F S Cc(G) should look something like 

^™ ; //I \ / / C(a^)C(y)-^(-2;+y)d2;dy = lim / / i{x) [ ^(y)(T^F)(y) dy dx , 

where the right hand side arises by using the usual trick from van Hove sequences and the 
compactness of the support of F. (That is, the difference of the two sides of the equation is 
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due to the difference between —An and —An + supp(F), which by the van Hove assumption 
is irrelevant in the Umit.) Of course in our case ^ is not a function of x and the integrands do 
not make sense. But the inner integral on the right-hand side is what should be N{TxF){S^) = 
N(F){T-x(,) and we can rewrite this 'autocorrelation' as 



n->oo lQ[An) Ja„ 



' An 

The term {T-xC){^) has no meaning. But Palm theory tells us how to go around this. We 



instead average over a small neighbourhood B of 0, and this brings us to Definition 3.13 



Theorem 3.14. Let G be a locally compact group whose topology has a countable basis and let 
{N, X, fi, T) be an ergodic spatial stationary process on G with second moment. The two-point 
correlation of ^ & X exists ^-almost surely and its value at F & Cc{G) is j{F). 

Proof. The function {N{1b)N{F)){T-xO is measurable as a function of x £ G and the 
Birkhoff theorem says that almost surely 

(9) hm — ^ [ {N{lB)N{F)){T.xOdx= [ N{1b)N{F) df, , 

n^oo lG{An) J a„ Jx 

meaning that the limit will exist and equal the right hand side. We shall prove that 



(10) 



lim 

_B-s>0 



N{lB)NiF) d/i = 7(F) 



X 



proving that the definition of Definition 3.13 works almost surely in ^ for each F G Cc{G) 
and that 7 does have the meaning of an autocorrelation. 

This has to be made to work simultaneously for all F in Cc(G), which will be shown in the 
usual way from a countable basis of Cc{G). This will prove Theorem 3.141 



Here are the details: It is enough to prove (10) for real- valued F. Using ([7]) we have for all 



Let F E 



F*Gd7 = {N{F) I N{G)) . 
and choose e > 0. Let i? be a measurable subset of G with compact closure. 



Then by Prop. 3.3, N{B) := N{1b) is defined and is a measurable L -function on X. By 



linearization we have 

/ N{lB)N{F)d^i 
Jx 



\b * F d'y 

lB{x)F{y - x) dx d7(y) 



Ib{x)F{x - y) dx d'y{y) 



Since F is uniformly continuous on G, for any sufficiently small neighbourhood B of 0, \F{x- 
y) — F{—y)\ < € for all x G .B and for all y £ G. Then 

I lB{x)F{x-y)dx-lG{B)F{-y) < [ 1b{x)\F{x - y) - F{-y)\dx < el^iB) 
JG JG 
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for all y G G, SO 



N{B)N{F)di^-k{B) / F(-y)d7(y) 



lBix)Fix-y)dxdj{y)-lGiB) / F(-y)d7(y) 



< 



< 



/ [ lB{x)\F{x-y)-F{-y)\dxd^{y) 



IrJB 



-elG{BMi-supp{F) + B) 



Since 7(y) = 7(— y) (7 is positve-definite and real, Prop. 3.1), we see that as S — (and 
correspondingly e — t- 0) we have a proof of (10). Together with ^ we have the desired 
interpretation of Def. 3.13 as the autocorrelation at ^ and its almost sure equality with j{F), 
and Prop. 3.14 is proved. □ 



Remark 3.15. This type of result was first established for certain uniformly discrete point 
processes on Euclidean space by Hof [19] (see [16] for the case of general locally compact 
abelian groups). This was then extended to rather general point processes by Gouere [T7] 
and to certain measure processes by Baake / Lenz [3]. A unified treatment was then given by 
Lenz / Strungaru [31]. Our result contains all these results (provided the underlying process 
is real). 

3.4. A second glance at second moments. The discussion above has shown that existence 
of second moments has strong consequences. It implies existence of the diffraction to dynamics 



map by Proposition 3.8 and further continuity properties by Corollary |3.2[ It turns out that 
a converse of sorts holds. This is investigated in this section. Along the way, we will also 
show that continuity of N implies an intertwining property of the map N and this will be 
crucial for our considerations. 



Any process {N, %, T) gives rise to two representations of C( 
The representation L lives on L^(G) and acts by 

LgF = G*F 



(and in fact even of L^(G)): 



for F,G € Cc(G). The representation T (extending the action T and therefore denoted by 
the same letter) acts by 



Tcf := I G{t)Ttf dUt) 



for G G Cc(G) and f € T-L. Continuity of N now yields an intertwining property: 
Lemma 3.16. Let {N,T-L,T) be a process with a continuous N . Then, 

NoLg = TgoN 

for allGe Cc(G). 
Proof. We have 

N{Lg{F)) = N{G*F) = N (^j^Git)TtFdlG{t) 
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By the linearity and continuity of A^, it commutes with taking integrals and 

n( [ G{t)TtFdlG{t)] = [ G{t)N{TtF)dlG{t)= [ G{t)TtN{F)dlG{t) 
\Jg / Jg Jg 

holds, which finishes the proof. □ 

We will need a somewhat stronger continuity property of A^. 

Definition 3.17. Let 1 < p, q < oo with 1/p + 1/q = 1 be given. The process {N,T-L,T) is 
said to be weakly (p, g)-continuous if for any compact K (ZG there exists a Ck with 

\{N{F)\N{G))\<Ck\\F\\ 



for all F,G G Cc{G) with support contained in K (see Cor. 3.2). 



Remark 3.18. (a) Note that continuity of is exactly weak (2, 2)-continuity. 

(b) By standard interpolation theory, we can conclude that weak (2, 2) continuity together 
with weak (1, oo) continuity implies weak (p, q) continuity for all p, q with 1 < p, (? < oo and 
1/p + 1/q = 1. 

Theorem 3.19. Let {N,T-L,T) be a stochastic process. Then, the following assertions are 
equivalent: 

(i) N has a second moment. 

(ii) N is weakly (p, q) -continuous for all 1 < p,q < oo with 1/p + 1/q = 1. 
In this case N is continuous and has the intertwining property that 

NoLg{F) =TgN{F) 

for aUF,Ge Cc(G). 



Proof. (i)^=^ (ii): This follows immediately from Corollary 3.2 

(ii)^=^ (i): As A^ is weakly (2, 2) continuous, it is continuous and the intertwining property 
follows from the previous lemma. In the remaining part of the proof we will only consider 



real valued functions. By Proposition |3.1[ it suffices to show existence of a measure 7 on G 
with 

{N{F) I N{G)) = 7(F * G) 
for all F,G €z Cc(G). A short calculation shows that for such a 7 the equality 

{N{F) I NiG)) = I F{y)ij*G){y)dh{y) 

must hold for all F,G £ Cc(G). The idea is now to 'reverse' this reasoning to conclude 
existence of 7. The main issue is to show continuity of the object Hq replacing the not yet 
defined 7 * G. This is shown using the intertwining property. Here are the details: We choose 
an arbitrary compact in G and assume without loss of generality that £ K. Let U be 
an open neighborhood of with compact closure and set Ki to be the closure of K + U. Set 
L = Ki-Ki + Ki -Ki. 

As A'^ is weak (l,oo) continuous, and L°°{L) is the dual of L^{L), we can find to each 
G € Cc(G) with support in Ki a function Hq in L^{L) with 

{N{F) I N{G)) = [ F{y)HG{y)dlG{y) 
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for all F G Cc(G) with support in L and 

(11) II^^gIIoo < C||G||oo, 

where C = C{L). As has the intertwining property, a short argument shows that 

Hg^*G2{x) = j Gi{y)HG2{x - y)dlG{y) 

for X G for all G2 S Cc(G) with support contained in K and all Gi G Cc(G) with support 
contained in U . This gives in particular, that Hgx*G2 is a continuous function (on K) for all 
such Gi,G2. Now (fTTl) yields 



\\{Hd*g — HG)\\cyo = \\Hd*g-g\\oo < C\\D * G — G||oo- 

Taking an approximate unit for D, we wee that Hq can be approximated uniformly by 
continuous functions. This shows continuity of Hg- Moreover, by (11) again the map 

G ^ Hg{x) 

is continuous for each x. Thus, we can indeed define the measure 7 with 

j*G{x) = Hg{x). 

By construction 7 has the desired properties. □ 



4. Almost periodicity 

Almost periodicity is closely linked to nature of diffraction via the Fourier transform. This 
concept allows one to compare properties of the autocorrelation and the diffraction measure. 
In particular, it can be used to characterize pure point and continuous diffraction. The 
considerations of this section play a role in subsequent parts of the paper: They are used 
in Section [5] to decompose an arbitrary processes into a part with pure point diffraction and 
a part with continuous diffraction, and in Section [6] to prove the equivalence of pure point 
diffraction spectrum and pure point spectrum. The material of this section derives from [15j 
and, in a more accessible account [50J, and from [301 [31] as well. For further studies of aspects 
of almost periodicity in our context we refer to [T71 HHl [23]. 

Almost periodicity is most commonly defined by a compactness condition: / G Cu(G) is 
almost periodic if the translation orbit G./ of / is compact. Here Cu(G) is the space of 
uniformly continuous C-valued functions on G. The key point, however, is in which topology 
is the translation orbit to be compact? In the case of the sup-norm topology the concept 
defines strong almost periodicity. In the case that the topology is defined by the family of 
semi- norms induced by the set of all continuous linear functionals on Cu(G), it is called weak 
almost periodicity. Strong almost periodicity coincides with H. Bohr's original concept of 
almost periodicity: / G Cu(G) is strongly periodic if and only if for every e > the set of 
t G Gforwhich||rt/-/|| sup < e is relatively dense. It is harder to get an intuitive feel for weak 
almost periodicity. Fortunately what one needs to know about it is fairly straightforward: 

• all positive definite functions on G are weakly almost periodic; 

• for every weakly almost periodic function / the mean of / defined as 

lim —I— I f{x + t)dlG{t) 
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exists for any van Hove (more generally F0lner) sequence {An} in G and any x G G 
and its value is independent of both x and the choice of van Hove sequence. 
One says that a weakly almost periodic function / is null weakly almost periodic if its mean 
is 0. These concepts lift to measures in a simple way: a Borel measure on G is called 
strongly (resp. weakly, null weakly) almost periodic if / * is a strongly (resp. weakly, null 
weakly) almost periodic function for every / G Cc(G). 

Positive definite measures on G turn out to be weakly almost periodic and are also, as we 
have noted already, Fourier transformable. 

Theorem 4.1. (Gil de Lamadrid, Argabright [15j) Every weakly almost periodic measure (j) is 
uniquely expressible as the sum of a strongly almost periodic measure and null weakly almost 
periodic measure. If the measure (p is Fourier transformable then so too are the strong and 
null weak components, and furthermore the Fourier transforms of these components are the 
pure point and continuous parts of (j). □ 

These considerations can be applied to strongly continuous unitary representations and 
hence to processes as well. This is discussed next. The crucial connection is given by the 
following (well known) lemma. 

Lemma 4.2. Let T be a strongly continuous unitary representation o/G onT-L then for any 
f £ H, the function i i— )• (/ | Ttf) =: Ff{t) is positive definite and hence weakly almost 
periodic and Fourier transformable. Its Fourier transform is the spectral measure pj-. 

Proof. As T is unitary, we have for any ti, . . . ,tn G G and ci, . . . , G C 

n n 

CkCiFf{tk-ti) = ||X]c«i^/(t/)f >0 

k,l=l 1=1 

and Ff is shown to be positive definite. The remaining statements follow from the above 
discussion. □ 

Combining the previous lemma and the previous theorem, we infer the following corollary. 

Corollary 4.3. Let T be a strongly continuous unitary representation of G on Ti then for 
any f £ Ti, the function t i— )• (/ | Ttf) is the sum of a strongly almost periodic function and 
a null weakly almost periodic function. This strongly almost periodic function is given by the 
(inverse) Fourier transform of the pure point part of the spectral measure pj and the null 
weakly almost periodic function is given by the (inverse) Fourier transform of the continuous 
part of the spectral measure pj . 

Of particular relevance is the question whether the spectral measures are pure point. Here, 
the following holds as shown in |31j . Lemma 2.1. 

Lemma 4.4. Let T be a strongly continuous unitary representation of G onT-L. Then, the 
following assertions are equivalent for f £%: 

(i) The map G — > %, t ^ Ttf , is almost periodic in the sense that for any e > the set 
{t G G : \\Ttf — /II < e} is relatively dense in G. 

(ii) The hull {Ttf : t G G} is relatively compact. 

(iii) The map G — > C, 1 1— )> (/ | Ttf), is strongly almost periodic. 

(iv) pf is a pure point measure. 

(v) / belongs to Tip. 
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The preceding equivalence allows one to show some stability properties of T-Lp. This is 
discussed next (see [31] as well for a proof of (a) of the following theorem). There, we say 
that a measure ^ is supported on a measurable set 5* if ^ of the complement of S is zero. 

Theorem 4.5. Let T be a strongly continuous unitary representation ofG on %. 

(a) Let C : % — > % be continuous with TfCf = CTff for each t € G and f & H. Then, 
C maps Tip into Hp. If f belongs to Tip and pf is supported on the subgroup S ofG, then so 
is pcf. 

(b) Let M : Ti X Ti — > % be continuous with TtM{f,g) = M(Ttf,Ttg). Then, M maps 
Hp X T-Lp into Tip. 

Proof The statements on C and M are clear from the equivalence of (v) and (ii) in the 
Lemma [4.4| and the continuity and intertwining property of C and M. The statement on the 
support in (a) is proven in jS^. □ 



We now further restrict attention to processes with second moment. Then, Theorem |4.1| 
can directly be applied to give the following result (which is well known from [5l I3 H [36 l I17j). 

Proposition 4.6. A stationary process {N, T-L, T) with second moment is pure point diffractive 
if and only if its autocorrelation is strongly almost periodic. It has continuous diffraction if 
and only if the autocorrelation is null weakly almost periodic. □ 

We can then combine the above results to obtain the following. 

Proposition 4.7. For any stationary process {N,T-L,T) with second moment the mapping 
t I—)- {N(F) I TtN{F)) is a weakly almost periodic function on G for all F G Cc{G). The 
process {N, Ti, T) then has pure point (resp. continuous) diffraction if and only if t ^ {N[F) \ 
TtN(F)) is a strongly (resp. null weakly) almost periodic function on G for all F G Cc(G). 
Furthermore, the diffraction is pure point if and only if for all F G Cc(G), t i— )• TtN(F) is a 
strongly almost periodic function on G with respect to the norm on T-L . 



Proof. The first statement follows from Lemma 4.2 We now turn to the second statement. 



By Theorem 4.1 and the definition of the diffraction spectrum we have pure point (resp. 



continuous) diffr actio n if and only if 7 is strongly (resp null weakly) almost periodic. Now, 



from Proposition 3.1 and a straightforward calculation we obtain for all F,G G Cjf(G) and 
ah t G G, 

(12) G*F* 7(t) = 7(G * T\F) = {N{G) \ TtN{F)) . 

With G set equal to F we then get the required almost periodicity properties of f 1— )• {N{F) \ 
TtN(F)) from the corresponding almost periodicity properties of 7, by our definition of these 
concepts. In the reverse direction, upon linearizing the expressions {N{F) \ TtN{F)) we 
obtain ( |12[ ). From the almost periodicity properties of 1 1— )• {N{G) \ TtN{F)) we then obtain 
the desired almost periodicity properties for all convolutions of the form G * F * 7(— t). An 
approximate unit argument allows us to get the strong (resp. null weak) almost periodicity 
of 7 from this, see [50], Corollary 9. Here, we use that these almost periodicity properties are 
stable under uniform convergence. 

The last statement is a consequence of Lemma |4.4[ This finishes the proof. □ 

We finish this section with a discussion of stability properties of almost periodicity in the 
case of spatial processes. We will need the following proposition on continuity of composition 
with a fixed function. The proposition is not hard to prove and can be found in [31j. 
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Proposition 4.8. Let {h„} be a sequence in L'^{X,fi) converging to h £ L^{X,fi). Let <f) 
he a continuous hounded function from the complex numhers to the complex numhers. Then, 
(j) o hn converges to (po h. 



Now, we turn to the fohowing consequence of Theorem 4.5 (see [3T] as weU). 



Corollary 4.9. Let {N,T-l,T) be a spatial process with % = L'^{X,fj.). Let i?i>, -0 : C — > C be 
bounded continuous functions and f,g£ Tip be given. Then, the following holds: 

(a) The function (f) o f belongs to %p and the spectral measure p^^^ is supported in the 

subgroup S of G if pj is supported in this subgroup. 

(b) The function {4>o f){ipog) belongs to Tip. If both pj and Pg are supported in the subgroup 

SofG then so is P^^of)i^og)- 

Proof. By the previous proposition the maps C : H — > %, f ^ (j) o f and M : % xTi. — > 
Ti, M(f, g) := ((/> o f){'ip o g) are continuous. They obviously commute with the action of G. 



Both (a) and the first statement of (b) fohow from Theorem 4.5 It remains to show the 
second statement of (b): Let be the subspace of %p consisting of elements whose spectral 
measure is supported on S. Then, it is not hard to see that T-Lp is a closed subspace of T-L. By 
(a) and the assumption on / and g, the spectral measures of both (j)o f and (f)og belong to Hp. 
Moreover, (f)o f and (j)og are bounded functions. Thus, it suffices to show that 'Hpr\L°°{X, p) 
is an algebra. This can be shown by mimicking the proof of Lemma 1 in [3] (see [25] as well). 
For the convenience of the reader we sketch a proof: 

Note first that every eigenfunction can be approximated by bounded eigenfunctions via a 
simple cut-off procedure viz if e is an eigenfunction, then for an arbitrary m > 0, the function 

(13) e(-)(x) := 1^^^)' 

I , otherwise 

is again an eigenfunction (to the same eigenvalue). Evidently, the e^"^^ converge to e in L^ 
as m — >• oo. Now, let non-zero functions a,b £ Tip n L°°{X,p) be given and choose e > 
arbitrarily. As a belongs to Hp, there exists a finite linear combination a = ^aiCi of 
eigenfunctions to eigenvalues in S with 



'^aiei\\2 < 



e 



By the conclusion after Eq. (13), we can assume that all are bounded functions. Thus, in 
particular, \\a ||oo < oo. 

Similarly, we can choose another finite linear combination h = '^bjCj of bounded eigen- 
functions to eigenvalues in 5 with 

lloo 

Then, 

||a6-a6||2 < ||a ||oo ||o - o lb + ||&||oo ||a - a ||2 < 2e. 



The product of bounded eigenfunctions to eigenvalues in S is again a bounded eigenfunction 
to an eigenvalue in S (as 5 is a group) and the claim follows. □ 
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5. Decomposing processes, the pure point-continuous split, and the Gelfand 

construction 

In this section we first discuss how any stationary process spUts into two subprocesses, 
one of which has pure point spectrum and one of which has continuous spectrum. If the 
original process has a second moment then so have the two subprocesses and their diffraction 



measure will be pure point and purely continuous respectively (Theorem 5.1). If furthermore 



the original process is spatial, then its pure point part is canonically isomorphic to a spatial 



pure point process (Theorem 5.6). This last result is based on the Gelfand construction of 
commutative Banach algebras. 

Whenever T is a strongly continuous unitary representation of G on the Hilbert space Ti 
and [/ is a T invariant closed subspace of Ti with corresponding orthogonal projection Pjj 
we can form the restriction Tu of T to U and this will be a strongly continuous unitary 
representation as well. The orthogonal complement C/"*" oi U in'H is also T invariant and we 
can form the restriction Tjj± as well. We will now focus on the special decomposition into 
the continuous and its point part. More precisely, we can define the continuous Hilbert space 
Tic and the point Hilbert space Tip by 

Tip := {f £ H : pf is a pure point measure}. Tic '■= {f £ Ti : Pf is continuous}. 

Then, both Hp and Tic are closed T-invariant with 

Ti = Tip © Tic- 

With the projections Pc and Pp onto Tip and Tic we then have Pc (B Pp = id and P^Tt = 
TtP^ =: T* for * = p, c. This gives the decomposition of the representation T asT = Tp@Tc 
with := P^:T, * = p,c. Obviously, all spectral measures associated to Tp are purely discrete 
and all spectral measures associated to Tc are purely continuous. Accordingly, we call Tp the 
pure point part of the representation T and Tc the continuous part of the representation T. 

We now turn to the case that we are given a stationary process {N, Ti, T). Of course, T is 
then a strongly continuous unitary representation of G on Ti. So, any T invariant subspace 
U gives rise to a decomposition of T. Moreover, the G-invariance of N then implies that A'^ 
can as well be decomposed into Nfj := PjjN and Njji_ = Pjj^N . Accordingly, the process M 
can be decomposed in to the two processes Njj := {Njj, U, Tu) and Afi/± := {Njj±, U-^,Tu±) 
and we write J\f = Mu ® J^u^ ■ 

As just discussed, T can be decomposed into its pure point part and its continuous part, 
into N = Np Q Nc with A'* = P*A, * = p,c. Hence, any process M := {N,Ti,T) can be 
decomposed into the two processes Mp := {Np,Tip,Tp) and Mc ■= {Nc,Tic,Tc) which are 
called the pure point part and the continuous part of M respectively. 

Now, assume that the original process has a second moment and associated autocorrelation 
measure 7 and diffraction measure = 7. By Theorem |4.1[ this 7 can be decomposed uniquely 
into its strongly almost periodic part "jsap and its null weakly almost periodic part ^owao and 
the Fourier transform of 7^0^ is the pure point part ujp of lo and the Fourier transform of 
"JOwap is the continuous part oJc of oj. It turns out that these are just the autocorrelation and 
diffraction measures of the pure point part and the continuous part of A^. This is the content 
of the next theorem. 

Theorem 5.1. Let M = {N,Ti,T) be a stationary process with a second moment and associ- 
ated autocorrelation 7, diffraction measure co and diffraction to dynamics map 6. Let ^sap o-nd 
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lOwap be the strongly almost periodic part and the null weakly almost periodic part ofj and ujp 
and ujc the respective Fourier transforms. Then, the point partAfp of M has a second moment 
with autocorrelation given by jsap o,nd diffraction measure given by Up and the continuous 
part Afc of M with autocorrelation given by 7o«,ap and diffraction measure given by uJc ■ 

Proof. By Proposition |3.1| it suffices to show that 

(14) 7,,p(F * G) = {PpN{F) I PpN{G)) and 7o«,ap(i^ * G) = {PcN{F) \ PcN{G)) 

for ah F,G G Cc(G). By linearity it suffices to consider the case F = G. Let Ft be the 
function TtF = F(- — t) and consider the function t i— )• ^{F * Ft). We will decompose this 
function in a strongly almost periodic part and a null weakly almost periodic part in two ways 
and then conclude the desired statement from the uniqueness of the decomposition. Here are 
the details: By definition of 'jsap and jowap we can decompose this function as 

7(F * Ft) = Jsap{F * Ft) + JOwap{F * Ft) 

with a strongly almost periodic function t i— >• ^sap{F * Ft) and a null weakly almost periodic 
function 1 1— )• jowapiF * Ft). On the other hand by definition of 7 and the fact that Pc and Pp 
are orthogonal we also obtain the decomposition 

7(F * Ft) = {N{F) I TtN{F)) = {PpN{F) \ TtPpN{F)) + {PcN{F) \ TtPcN{F)). 

Now, t I— )• {PpN{F) I TtPpN(F)) = J{t,$,)dpp 7v(f) is strongly almost periodic as it is 
the inverse Fourier transform of a pure point measure and t 1— )• {PcN{F) \ TtPcN(F)) = 
J {t,C)dpp^N{F) is null weakly almost periodic as it is the inverse Fourier transform of a con- 
tinuous measure. From Theorem 14. II we conclude that 

lsap{F * Ft) = {PpN{F) I TtPpN{F)) and 7o»ap(F * Ft) = {PcN{F) \ TtPcN{F)) 

for all t E G. Setting t = we obtain □ 

In the situation of the previous theorem we can easily see that the diffraction to dynamics 
map 6 : L^(G,a;) — > V. can be decomposed as 

(15) e = ep® Oc. 

Here, Op and 9c are the diffraction to dynamics maps of Np an Nc respectively and and the 
natural decompositions L^(G,a;) = L^(G,a;p) © L^(G,a;c) and T-L = Tip ® TLc are implicitly 
used. 

We now restrict attention to spatial processes. In this situation we can induce subprocesses 
via subalgebras of L°° . This is discussed next. We will use the main result of Gelfand theory 
that any commutative C*-algebra A with a unit is canonically isomorphic to the algebra 
C{X_a) of continuous functions on the compact space given by the maximal ideals of A. 
We will denote the space of all continuous bounded complex valued functions on the complex 
plane by Cb(C). 

Proposition 5.2. Let N := {N,Ti,T) be a spatial process with Hilbert space Ti = L^(X, /i). 
Let A be a G-invariant subalgebra of L°°{X,fj,) which is closed under complex conjugation, 
contains the constant functions and is closed with respect to the sup norm. Let U = U{A) 
be the subspace of Ti generated by A and let X_4 the maximal ideal space of A. Then, Xj\^ 
can be equipped in a unique way with a G action and a G-invariant measure p^ such that the 
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canonical Gelfand isomorphism J : C(X_4) — > A extends to an unitary G-map (also called 
J) 

J:L^{Xa,I^a)^U. 
The map J is compatible with the algebraic structure in that it satisfies 

• Jifg) = J{f)J{9) for all f G L~(X^,^^) and g G L'^{Xj^,hX) and 

• J{(p{f)) = (j){J{f)) for all f ^ LP' and for all continuous bounded (j) : C — > C. 

The subspace U is a G-invariant subspace and N can be decomposed as 

where Mu is isomorphic via J to the spatial process Na with Na = J-^PuN, Ua = 
L'^{Xa, fJ-A)- If "is ergodic, so isMu- 

Proof. All statements are rather straightforward up to the compatibility of J with the alge- 
braic structure and the ergodicity. 

This compatibility with algebraic structure will be discussed next: By definition, J is 
an algebra isomorphism from C{Xa) and A. In particular, J{fg) = J{f)J{g) for f,g G 
C{Xa)- By a simple approximation argument, we then obtain J{fg) = J if) Jig) for all 
/ G L°°iXA) and g G L^(X_4). Here, we use that {fgn} converges to {fg} in L'^HXa) 
whenever {gn} converges to g in L^(X_4) and / is bounded. This shows that J is compatible 
with multiplication. 

As J is an isomorphism of algebras, we also have J{qif)) = qiJif)) for any polynomial q 
and any / G C{Xa)- By a Stone/ Weierstrass type argument, this gives J((/)(/)) = (/>(J(/)) 
for any (p G Cb(C) and / G C(X_4). Consider now the set 

C := {/ G L\Xa,Ixa) ■■ 4>iJ{f)) = JiHf)) for all ^ G ^(C)}. 

Then, C(X_4) C C by what we have just shown. Moreover, C is closed in by Proposition 
4.8 Thus, C must agree with L^(X_4,//_4) as C{Xa) is dense in L^(X_4,//_4) (as {(pihn)} 



converges to (pih) whenever {/in} converges to h). 

It remains to show the statement on ergodicity. By definition ergodicity means that the 
eigenspace of T to the eigenvalues 1 is one-dimensional. This is obviously stable under the 
above constructions. □ 

Definition 5.3. The process Ma constructed in the previous proposition is called the sub- 
process of M induced by A. 

Theorem 5.4. Let M := (A^, "HjT) be a spatial process with Hilbert space T-L = L^(X, //). Let 
A be the closure with respect to the supremum norm of the subalgebra ofL°°iX,fi) generated 
by ipi o A^(Fi) . . .Tpn o N{Fn) with n G N, Fj G Cc(G) and il^j continuous bounded functions 
from C to C. Let U = UiA) be the subspace ofH generated by A. Then U contains the range 
of N and the subprocess Mu is isomorphic to the full spatial process Ma- The process Ma is 
ergodic if M is ergodic. Moreover, 

(NAiF) I NAiG)) = {NiF) \ NiG)) 

holds for all F,G £ Cc(G). In particular, M has a second moment if and only if Ma has a 
second moment and their autocorrelations agree in this case. 
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Proof. By definition of the range of N is obviously contained in the subspace U . Thus, 
PjjN = N and therefore 

{Nu{F) I Nu{G)) = {N{F) \ N{G)) 

holds for all F,G € Cc(G). Given this all statements of the theorem follow immediately from 
the previous proposition. □ 

Remark 5.5. The theorem says that to any (ergodic) stationary process M we can find an 
(ergodic) stationary subprocess Mjj which is full and has the same second moment. In the 
situation of the theorem fullness of the original process J\f then just means that J\f = Mjj- 

Now let {NjTi, T) be a spatial process with H = L'^{X, ^). Let V C L'^{X^ fi) be the closure 
of the linear span of products of the form 

V-i o Np{Fi) ...i^no NpiFn), n e N, € ^(C), € Ce(G). 

Then, y is a closed invariant subspace of L'^{X, fi). Let Py be the orthogonal projection onto 
V. Then, Py commutes with T. Then, {Ny, V, Ty) with Ny = PyN and Ty = PyT = TPy is 
a stationary process. We call {Ny, V, Ty) the augmented point part of the stationary process 
(Nj'HjT). Let V-^ be the orthogonal complement of V in Ti. Then, V-^ gives rise to a 
stationary process ( A^y± , y-*- , Ty ± ) and N is the sum of these processes. 

Theorem 5.6. Let M := {N,T-l,T) be a spatial process with Hilbert space Ti = L'^{X,fj.) and 
associated augmented point part My := {Ny, V, Ty). Then the following holds: 

• M = N'y® Myi^ ; 

• My is isomorphic to a full spatial process; 

. {PyN{F) I PyN{G)) = {PpN{F) \ PpN{G)) and{Py^N{F) \ Py±N{G)) = {PcN{F) 
PcN{G)) hold for all F,G e Cc{G). 
If M has a second moment with autocorrelation 7 and diffraction u furthermore the following 
is true: 

• My has a second moment with autocorrelation ^sap o-nd diffraction measure given by 
the pure point part of u. 

• My± has a second moment with autocorrelation jowap o,nd the diffraction measure 
given by the continuous part ofu. 

Proof. Let Mp = {Np,'H,T). Let the algebra A be the closure with respect to the sup-norm 
of the linear span of products of the form 

V-i o Np{Fi) . . . o Np{Fn), n e N, € ^(C), € G,{G). 

Then, V is just the subspace generated by A. 

We can then apply the previous theorem to Mp and the algebra A. This gives a decompo- 
sition Mp = My © Myi^ , where My is isomorphic to a full spatial process and 

{Ny{F),Ny{G)) = {Np{F),Np{G)) = {PpN{F) \ PpN{G)) 
holds for all F,G G Cc(G). Here, the last equality follows by definition of Np as PpN. 

It remains to show that PyN{F) = PpN{F) (then all remaining statements follow easily). 
To do so, we first note that the range of Np = PpN is contained in the range of U by 
construction. Moreover, by general principles the algebra A and hence the subspace V is 
contained in the pure point part of T [31 [25]. Thus, PyN = PpN holds. 

Similar remarks apply to iVc- Note that Py±N = N - PyN = N - PpN = P^N . □ 
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The considerations of this section and in particular the preceding theorem have the following 
consequence: whenever we are given a spatial process we can restrict attention to its point 
part and obtain a full spatial process with pure point spectrum. In particular, the whole 
theory developed below for full spatial processes with pure point spectrm will apply to the 
point part of any spatial process. 

6. Spatial processes with pure point diffraction 

In this section we discuss spatial processes with pure point diffraction. It is these processes 
that we will be able to classify. We show that for these processes the two notions of pure 
pointedness of the spectrum defined above viz pure point diffraction spectrum and pure point 



dynamical spectrum actually agree. As noted in Theorem 5.6, any spatial process with second 
moment can be decomposed into a full spatial process with pure point diffraction and a general 
stationary process with continuous diffrction. Thus, the results of this section apply to the 
corresponding parts of any spatial process. 

Let (A^, X, /X, T) be a stationary process with pure point diffraction as discussed in Section 
[2| Thus, its diffraction measure cj is a pure point measure and the set of its atoms denoted 
by 5 = 5(a;) is given by 

5 = {A; G G : u){h) > 0}. 

Let 

e : L^{G,uj) — > L^{X,m) 
be the associated diffraction to dynamics map and := 9 {Ik), k £ S. 

While the subspace 9{L'^{G,uj) of L'^{X,fi) may be small in some sense, in another sense 
it controls the whole Hilbert space {X, n) due to our assumption of pure point spectrum. 
A precise version is given next. 

Theorem 6.1. Let M = {N,X, fi,T) be a spatial stationary process, which is full and pos- 
sesses a second moment. Then, the following assertions are equivalent: 

(i) The diffraction measure uj of (N, X, n, T) is a pure point measure. 

(ii) The representation T has pure point spectrum. 

In this case, the group £ of eigenvalues ofT is generated by the set 5(a;), and any eigenfunc- 
tion is a (multiple of a) product of eigenf unctions of the form fk, k ^ S{uj). 

Proof. We mimick the argument of |31] (see [3l [25] as well). By definition of the spectral 
measures pf given in ^ we have 



{N{F) I TtN{F)) = / (7,i)rfP7V(F)(7) 
Jg 

for all F G Cc(G). Moreover, by definition of uj and 7 and the fact that A is a G-mapping 
we find 

(A(F) I TtN{F)) =jiF*FXT^f )) = f{-f,t)\F\^duij) 



first for real- valued F £ Cc(G) and then by linearity for all F S Cc(G). Taking inverse Fourier 
transforms we infer 

(16) pNt^F) = \Ffduj 

for all F G Cc(G). 
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The implication (ii) =^ (i) is now clear from (16). (We also saw this earlier from the 
properties of 9.) We next show (i) =^ (ii): As the diffraction measure is a pure point 
measure, the spectral measures associated to N{F) are pure point measures (supported on 



S) by (16). Thus, by part (a) of Corollary|4.9[ the spectral measures of ^oA^(F) for ip £ Cc(C), 



F G Cc{G), are then also pure point measures supported on the group generated by S. By 



the fullness assumption and part (b) of Corollary 4.9 we then infer that T has pure point 
spectrum with eigenvalues contained in the group generated by S. 

Now, multiplying eigenfunctions fk, k € S, we obtain eigenfunctions for arbitrary elements 
in the group generated by S. This shows that the group generated by S is indeed the group 
of eigenvalues. At the same time it gives an eigenfunction for each eigenvalue. By ergodicity 
the multiplicity of each eigenvalue is one and these eigenfunctions are all eigenfunctions. □ 



Remark 6.2. (a) There is quite some history to Theorem 6.1: The implication (ii) =^ (i) 
is sometimes known as Dworkin argument. It was first established by Dworkin [13] with later 
extensions by Hof [19] and Schlottmann [46] (see [3Hll9] for remarkable applications as well). 
In special situations equivalence was then shown by Lee/Moody/Solomyak in [25j. This was 
then generalized in various directions by Gouere [T7], Baake/Lenz [3], Deng/Moody |10) . and 
Lenz/Strungaru |31j (the latter result containing all earlier ones). The above result contains 
all these results (provided the underlying process is real). 

(b) The equivalence between diffraction and dynamical spectrum cannot hold for the com- 
plete spectrum, as discussed by van Enter/Mieskicz [H]. Indeed, [14j gives an example of 
mixed spectrum where the pure point part of the dynamical spectrum is completely missing 
in the diffraction (up to the constant eigenfunction). 

Due to the previous theorem we do not need to distinguish between pure point diffraction 
and pure point dynamical spectrum when dealing with full processes with second moments. 
This suggests the following definition. 

Definition 6.3. A spatial stationary process (N, X, fi, T) is called pure point if it is full, 
possesses a second moment and its diffraction measure is a pure point measure. 



If G 5 in the situation of Prop. 3.12| then by Lemma 3.11| /o is a real function belonging 



to the k = 0-eigenspace of L'^{X,m) and ||/o|| = w(0)^/^. Since constant function 1 x '■ ^ ^ 
for all ^ G X is an eigenfunction for k = and the dynamical system is ergodic, it spans the 
k = 0-eigenspace of L^(X, m). Thus /o = ±ijj{0)^^'^ Ix- Now, we have already pointed out 
that if is a stationary process then so is any non-zero real multiple of it. In particular ibA^ 
are stationary processes and one of them has the corresponding function /o = uj(Oy^'^ Ix- It 
causes unnecessary awkwardness later on in the discussion of phase forms to deal with the 
case /o = — a;(0)^/^ix> so wish to normalize A^ to avoid this situation: 



Assumption 6.4. 



In the pure point ergodic case we shall always assume that 

fo = oj{oy/^ix. 



7. Relators and associated phase forms 



As mentioned already, our aim is to describe all stationary processes with a given pure 
point diffraction measure u. This description will be given in terms of a set of objects called 
relators arising out of the Bragg spectrum of uj. While the later sections need the material of 
this section, it can be read independently of the previous ones. 
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Let G be a locally compact abelian group. Let 5 C G be given with S = —S and let £ 
be the subgroup of G generated by S. We shall give this group the discrete topology, and in 
order to keep this clear we shall usually write £d instead of £. 

In the sequel we often use m-tuples {ki, . . . , km) £ 5™ (for various m) of elements of S or 
of £d and we often simply denote these by bold letters k. If k = {ki, . . . , km), I = {h, ■ ■ ■ ,ln) 
then kl is the concatenation 

kl = {kl, . . . , km, Wi • • ■ , In) 

and 

[k] := /ci H h G £d- 

In this way S := IJm=o'^™ becomes a monoid under concatenation, with the empty 0-tuple 
as the identity element. 

A relator is any m-tuple k = [ki, . . . , km) £ 5™" with 

[k] = ki + --- + km = 0. 
The empty tuple is also taken to be a relator. Thus the relators are a subset of S. Let 
Zn := {k=(A;i,...,fc„)G5" : [k]=0},n>l; 
Zo := {0}; 

oo 

Z := \JZn. 

n=0 

Note that Zi = {(0)} if £ S , otherwise it is empty. 

Concatenations of relators are relators and this makes Z a submonoid of S. 

We introduce an equivalence relation on S by transitive extension of the three rules: 

(Rl) k ~ 1 if 1 is a permutation of the symbols of k; 

(R2) if k = (fei, . . . , km) S S and € 5 then {ki, . . . , km,0) ~ (/ci, . . . , km) (along with the 
first rule, this means that when £ 5, zeros can be dropped or added wherever they 
appear) ; 

(R3) k ~ 1 if 1 can be obtained from k by inserting removing pairs {k, —k}, k € S. 
The second item here is related to our Assumption |6.4[ see below. 

It is easy to see that k ~ 1 , k' ~ 1' ^ kl ~ k'l', so multiplication descends from S to 
S~ := S/ ~, whereupon S~ becomes an Abelian group under this multiplication with 0~ as 
the identity element. Note that if G 5 then 0~ = 0~. We give S~ the discrete topology, so 
in particular it is a locally compact Abelian group. 

Since the sum [k] = ki + ■ ■ ■ + km of an element of 5*" is constant on the entire equivalence 
class k~, we see that Z is the union of all equivalence classes with component sum equal to 0, 
and we obtain Z := Zj ~ as a subgroup of S~. In addition, we see that there is a surjective 
homomorphism 

(/> : S~ — > £d, with k ^ [k] 

and its kernel is precisely Z. 

Definition 7.1. The group Z = Z{S) is called the relator group. The elements of its dual 
group Z are called phase forms. Thus a phase form is a group homomorphism 

a* :Z ^ U{1) 

(the latter being the unit circle in C). 
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Remark 7.2. (a) It is useful to think of relators in terms of the cycle structure of the Cayley 
graph of the group of eigenvalues with respect to the generators S. In this way the group Z 
could then be seen as a a kind of abelianized homotopy group or a first cohomology group of 
the graph. This works in the following way: We begin with £ and the set S which satisfies 
S = —S and generates £. The Cayley graph has vertices £ and edges consisting of all pairs 
(x, X + k) where x £ £ and k G S. Since for each edge (x, x + k) we have the reverse edge 
{x + k, X + k — k = x) , we can think of the graph as being non-directed. Since S generates £, 
the graph is connected. Each k = (fci, . . . , km) can be thought of as an edge sequence. Given 
any x £ £ we obtain a path x, x + ki,x + ki + k2, ■ ■ ■ , x + ki + ■ ■ ■ -\- km from it. This path is 
a cycle if and only if k is a relator, that is, [k] = 0. 

There is a sort of homotopy of edge sequences, which is generated by rules corresponding 
to (Rl), (R2), (R3) above. The first amounts to saying that for all x £ £ and for all k,l £ S, 
the paths x,x + k,x + k + l and x,x + l,x + l + k are homotopic; the second says that if G 5 
then the path x,x + is homotopic with the empty path from x; and the third says that the 
path x,x + k,x + k + (—k) is homotopic to the empty path from x. The group Z may be 
thought of as a sort of homotopy group for cycles originating (and terminating) at 0. 

(b) It should be noted that the homotopy group Z may be infinite even when the group 
G, and hence the associated graph, is finite. We shall see this in the examples below. 

The phase forms play a central role in homometry problem. We note that they depend 
only on 5, not on the actual values of the diffraction measure uj. We shall find it conve- 
nient to sometimes abuse the notation and write things like a*{ki, . . . ,km) where we mean 
a*{{ki, . . . , km)^) when using the phase form a* . 

Of the various k = {ki, . . . , km) that can represent a given element in Z there is (at least) 
one of minimal length n. This minimal length is denoted by len(k), and is called the reduced 
length of k. Here, of course, the length of {ki, . . . , km) is given by m. 

Define Zn to be the set of elements which have reduced length less than or equal to n. We 
have 

oo 

Z = [JZn, 

n=0 

ZnZp C Zn+p for all non-negative integers n and p. 

We let J^{S) be the Abelian group (with discrete topology) generated by symbols e{k), 
k £ S, subject only to the relation^ e(0) = 1 (if S 5) and e{—k)e{k) = 1 for all k £ S. 

Proposition 7.3. The mapping ip : e{k) i— t- (k)^ for all k £ S extends to an isomorphism 
if : T{S) S~. 

Proof. The existence of (/9 as a surjective homomorphism is clear since S~ is the Abelian 
group freely generated by the symbols (fc)~ subject only to the relations defining J-{S). 
These correspond precisely to {k)^{—k)^ = (k, —k)^ = 0~ and (0)~ = 0~ when E 5 which 
hold in S~. 

A typical element e(fci) . . . e{km) of -F(5) is in the kernel of ^ if and only if (A^i) . . . (^km) — 
{ki, . . . , fem)~ = 0- This happens only if the non-zero components of {ki, . . . , km) cancel in 
pairs k, —k. Then also the corresponding terms in e(fci) . . . e{km) cancel in pairs. Since also 
e(0) = 1 in T{S) when G 5, we see that e{ki) . . . e{km) reduces to the identity element. □ 

^The second of these relations comes from the underlying assumption that we are dealing with real stationary 
processes. If this condition were dropped then these relations would also be dropped. 
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We find it convenient to identify J-{S) and S~ through Prop. 7.3 Having done this we 
obtain from (p : S~ — > S^i, defined above, the homomorphism 

(p : T{S) — frf, with e{k) i— )• k , 

whose kernel can be identified with Z. 

Thus, we have the following exact sequence of groups (all given the discrete topology) 

(17) i^z^ J^{S) ^Sd^l. 
Dualization then gives the exact sequence 

(18) 1^ — Zi — f{S)i — — 1, 

where T := S^is & compact Abelian group by Pontryagin duality. Note that we get surjectivity 
from the continuity of the mappings and the compactness of these groups. In fact, surjectivity 
of the map T{S) — Z is crucial for our subsequent considerations. 

Definition 7.4. The mapping a : S — > U{1) satisfies the m-moment condition if 

(19) a{ki)a{k2)...a{km) = l 
whenever ki, . . . ,km ^ <S and ki + • • • + km = 0. 



We shall see the reason for calling these moment conditions in §15.1 Notice that the first 



moment condition is empty (and so trivially holds!) unless S 5, in which case it says that 
a(0) = 1. The second moment condition is equivalent to the statement that a{—k) = a{k) for 
all k G S. If € 5 then the second moment condition alone gives a(0)a(0) = 1, so a(0) = ±1. 



See Remark 6.4 If the first moment condition is not empty and holds then a(0) = 1 and then 
the m-moment condition includes the n-moment condition for all n < m since additional slots 
can be filled with zeros. 

Obviously, any mapping a : S — > U{1) satisfying the first and second moment conditions 

determines a character J-^{S), i.e. an element of J^{S), via e{k) i— )• a{k) for all k € S. 
Conversely every character of J-{S) clearly determines a mapping a satisfying the first and 
second moment conditions. (Notice that when G 5 then e(0) is the identity element of 
7^(5) and also a(0) = 1.) In the sequel we will use these two concepts interchangeably and 
use the symbol a both as a mapping on S and as the corresponding character on J^{S). When 

thought of as elements of -F(5) they are called elementary phase forms. 

By restriction any elementary phase form a determines a phase form a* , that is, an element 



of Z. The exact sequence (18) shows that every phase form arises by restriction from an 
elementary phase form and two elementary phase forms restrict to the same phase form if 
and only their ratio is in T, that is to say, if and only if their ratio is a character on £d- This 
is the essential part of determining equivalence of pure point processes. 

Proposition 7.5. (i) Any mapping a : S — > U{1) satisfying the first and second mo- 
ment conditions determines a unique elementary phase form and every elementary 
phase form arises in this way. 

(ii) Any phase form is the restriction of an elementary phase form. The ratio of any two 
elementary phase forms giving the same phase form is an element of T. 

(iii) An elementary phase form satisfies the m-moment conditions for m = 1, . . . ,n if and 
only it kills Zn- 
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(iv) An elementary phase form a determines a character on Sd; '^•e. lies in T, if and only 
if it satisfies all m-moment conditions, m = 1,2, ... . 

Proof. Parts (i) and (ii) are already done above. 

Zn is generated by all the k = {ki, . . . , km) with ki + ■ ■ ■ + km = 0, m < n. An elementary 
phase form a kills the equivalence class of k in J-^{S) iff 

a{ki) . . . a{km) = 1 • 

This proves (iii). An elementary phase form a kills all of Z iff a* is trivial which happens iff 
a = which gives (iv). □ 

Of particular interest to us is the case where the group Z is generated by Zn for some 
n > 2. This can be understood on the level of phase forms in the following way. 

Lemma 7.6. Let n > 2. Then the following assertions are equivalent: 

(i) {Zn) = Z. 

(ii) Any mapping a : S — > U{1) satisfying the m-moment conditions for 1 < m < n 
extends to a character on £d (i-s. lies in T). 

Proof. (ii)^=^ (i): The fact that a : S — > U{1) satisfies the m-moment conditions for 1 < m < 
n is equivalent to saying that the corresponding character a on J-{S) kills Zn, or equivalently, 



kills the subgroup (Zn) that it generates. But Prop. 7.5 says that a extends to a character on 
£d iff CL kills Z. Suppose that fact that the mapping a : S — > U{1) satisfies the m-moment 
conditions for 1 < m < n is enough to sufficient to guarantee that whenever a kills {Zn) it 
must also kill Z. Then {Zn) = Z because {Zn) C Z, both are closed subgroups of J-^{S), 
and the Pontryagin duality theory says that the characters can distinguish distinct closed 
subgroups. 

The implication (i)=^ (ii) is clear. □ 

The preceding considerations suggest to consider 

no := uq^Z) := inf{n £ N : Z generated by Zn}. 
Here, the infimum of the empty set is defined to be oo. 

Remark 7.7. If £ = S + ■ ■ ■ + S (with r summands) then any relator can be written as 
a product of relators of length at most 2r -|- 1 so Z2r+i generates Z. Thus, in this case 
no<2r + l [29\. 

Definition 7.8. The restriction of a phase form or an elementary phase form to Zm is called 
its mth- moment. 

Proposition 7.9. Two elementary phase forms a and b give rise to the same mth moments 
if and only if their ratio u = b/a satisfies the mth moment condition. 



Proof. Elementary phase forms a, b give rise to equal mth moments if and only if a{ki , . . . ,k„ 
b{ki, . . . , km) whenever ki, . . . ,km G 5 with ki + ■ ■ ■ + km = 0- However, being elementary 
phase forms, these equations can be written as a{ki) . . . a{km) = b{ki) . . . b{km), which leads 
to the equivalent statement that 

u{ki) . ..u{krn) = 1 

whenever ki, . . . ,km G 5 with ki + ■ ■ ■ + km = 0. This is the mth moment condition. □ 
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Two elementary phase forms a, b give rise to the same phase form a* if and only if their 



ratio is an element of T and so satisfies all the moment conditions. By Prop. 7.9 this means 
that all their moments are the same, so that there is no ambiguity in speaking of the mth 
moments of phase forms. 

Lemma 7.10. Assume (2^) = Z for some n > 2. Then two phase forms agree if and only 
if they have the same mth moments for m = 1, . . . ,n. In particular, two elementary phase 
forms restrict to the same phase form if and only if they have the same mth moments for 
m = 1, . . . ,n. 

Proof The first statement is clear. The second statement is a direct consequence of the first 
statement. □ 

Remark 7.11. Any of the phase forms a that we are interested in from the point of view 
of diffraction satisfy the first and second moment conditions by assumption. The relevance 
of the previous lemma is the following. If two stationary processes have the same pure point 
diffraction and their associated group of relators satisfies {Zn) = Z, then the processes are 



isomorphic if they have the same mth moments for m = 1, . . . ,n (see Cororllary 9.2 and 



Corollary 10.3 for precise statements). 



8. From spatial stationary processes with pure point diffraction to phase 

FACTORS 

In this section we assume that we are given an ergodic spatial stationary process with pure 
point diffraction. We will discuss how this process gives rise to a diffraction measure lo, a set 
S which generates a group £, canonical eigenfunctions fk, k £ S, and a phase form a* . In the 
next section we will see that (w, a*) in some sense uniquely determines the process. 

Let {N, X, n, T) be an ergodic stationary spatial process with pure point diffraction. From 
the considerations of Section [2] we then obtain the following: The diffraction measure uj of 
(A^, X, //, T) is a pure point measure and the set of its atoms, which are called Bragg peaks, 
is denoted by 5 = 5(a;) i.e. 

S = S{oj) := {A; G G : uj{k) > 0} 

and satisfies S = —S. We let £ be the subgroup of G generated by S and let £d denote 
this group when it is given the discrete topology. We are now clearly in the situation of the 
previous section. In particular, there is an associated relator group Z = Z{ijj) at our disposal. 
Let 

e : L^{G,uj) — > L^{X,m) 

be the associated diffraction to dynamics map and for each G G, let 1^ be the characteristic 
function on the set {A;}. Let U be the action of G acting on L^(G,cl!) defined by 

Uth{k) := {k,t)h{k) 

and observe that whenever Lj{k) 7^ 0, 1^ is a A; eigenfunction for this action,. It is easy to 
see that up to scaling these are the only eigenfunctions of L'^{G,uj). By the orthogonality of 
eigenfunctions we have, 

(Ifc I l/c')g = w(/c)(5fc,fc' 
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and 

{F\G)^ = Y,mGikMk) 
kes 

for all F,G £ L^{G,u)). 

To each element of 5 there exists a corresponding canonical eigenfunction. More specifically, 
we define fk, k £ S, hy fk = 0{lk). Then each fk is an eigenfunction in L^(X, ^u) for 
the eigenvalue k and ||/fc||2 = = uj{k)^^'^. Since \fk\ is a constant function (due to 

ergodicity) and is a probability measure, \fk\ = f^(A;)^/^. Then, for all ki, . . . ,km £ S, 

\\fki---fkj\^= / fk,---fk^fk,---fkm'^f^= / \fki\^ ■■■\fkj^dfi = uj{ki)...uj{km). 



At this point we invoke Assumption 6.4 which says that if G 5 then /o = oji^fl'^lx. 
For k = (A;i, . . . , km) with kj G S and [k] = the product fki ■ ■ ■ fk^ is an eigenfunction to 
and hence (by ergodicity) a multiple of Ix, which yields 

(20) fk,...fk^=a*{ki,...,kmMki)^/\..uj{kmf/^ 

for some a*(k) = a*{ki, . . . , km) G U{1). In this way we have a mapping a* : Z = Z(S) — > 
U{1). It clearly r espec ts the multiplication by concatenation on Z. In view of our assumption 



on /o and Lemma 3.11 we have f-k = fk, which shows this a* is well-defined on the equivalence 



classes of Z = Z{S) defined in ^ and in this way we obtain a phase form 

(21) a* : Z ^U{1). 

The preceding considerations show that any pure point stationary spatial process comes 
with a natural measure uj and a phase form a* as summarized in the following proposition. 

Proposition 8.1. Each ergodic pure point stationary spatial process (N, X, iJ,,T) gives rise 
to a pair (w, a*) consisting of a pure point measure u on G characterized by 

[ FGdu = {N{F) I N{G)) 
Jg 

for all F,G £ Cc(G) and the phase form a* on the relator group Z = Z(u) satisfying 

fk,... fk^ =a*{k,,..., km)0j{klY'^ . . . u{km)^/\ 

Next we show that u and a* are a complete set of invariants for the underlying process. 

9. Isomorphism of pure point processes 

Recall that we have introduced a notion of isomorphism between spatial processes in ^ In 
this section we show that two pure point ergodic spatial stationary processes are isomorphic 



(in the sense of Definition 2.5) if they yield the same diffraction measure uj and the same 
phase form a* . 

We remind the reader that our concept of a pure point spatial stationary process given in 
Definition |6.3| entails that the process in question is full and possesses a diffraction measure 
(which is pure point). 

Theorem 9.1. Two pure point ergodic stationary spatial processes based on the same group 
G are isomorphic if and only if their associated diffraction measures and phase forms are the 
same. 
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Proof. Assume that the full processes {N, X, fi,T) and {N' , X' , fi' ,T') are pure point ergodic 
and have the same diffraction measure lo and the same phase form a* . Set S := S{uj) = S{uj'). 
Let fk G L'^{X,fi) and G L'^{X' , /j,'), k £ S, he the corresponding natural eigenfunctions. 

Thus, N{F) = Y.kes^ik)fk and N'{F) = EkeS^Wf'k ^or all F G Cc(G). Moreover, by 
Theorem 6.1 products of the fk and the respectively provide orthonormal bases consisting 



of eigenfunctions of L^(X, /i) and L^(X',/i') respectively. 

We construct an isomorphism M as follows: We define M to be the unique linear map with 

fki ■ ■ ■ fkm ^ fki ■ ■ ■ fkm ■ 

This mapping is well-defined, for suppose that fk^--- fk^ and fi^... fi„ are linearly dependent. 
Then ki + ■■■ + km = h + ••• + In, so ki + ■■■ + — h — ••• — In = and 

Ai • • • fkmf-h ■ ■ ■ /-/n = ^(^l)^^^ • • • io{km)ujih)^^'^ . . . L0{lny^'^a*{ki, ...,km,-h,..., -In) ■ 

Using Lemma [3.11| we obtain 

fk,... fk^ = uj{kif'^ . . . oj{kmfl^u:{lir^/^ . . . u{Q''/^a*{ki, . . . ,km, -h, . . ■ , -ln)fw ■■■fir.- 

By our assumptions we obtain the same relation when the fs are replaced with /'s, so the 
same linear dependence occurs in 

As products of the fk and respectively give an orthogonal basis of L?'{X, ji) and Lp'{X' , fi'), 
and as 

\\fk,...fkjf = u;{k,)...u;{km) = \\f'k,...f'kjf. 
our mapping is a unitary map. In particular it is continuous. 

The definition of M easily shows that it intertwines T and T' and maps N{F) to N'(F) 
for all F £ Cc(G). Moreover, by definition M satisfies 

Mifk, . . . fkJh ■■■fiJ = Mifk, . . . fkjMifi, ...fij 

for all ki,... ,kn G 5 and h,. . .,lm G S. Note that both fk, ■ ■ ■ fk„ and M{fk, ■ ■ . fkj = 
f'ki ■ ■ ■ f'kn bounded functions. Thus, we can take linear combinations (in L?) and their 
limits to obtain 

M{fk,...fk^g) = M{fk,...fk^)M{g) 

for all g G L'^{X, fi) and ki . . .k^ £ S. In particular, 

M{fk,...fk^9)=M{fk,...fkjM{g) 

holds for all g G L°°(X, and ki . . .kn G S. Another approximation argument now yields 
M{fg) = M{f)M{g) for all f,g £ L°°{X,iJ.). This proves the isomorphism. 

Conversely suppose that M is an isomorphism between the pure point ergodic stationary 
spatial processes (A'', X, /x, T) and (A^', X' , /u', T') on G. As M is a unitary mapping, one sees 
from directly from the definitions that their correlation measures, and hence their diffraction 
measures are the same (if second moments exist), and then that the diffraction to dynamics 
mappings are related hy 6' = M o 8. Then looking at the eigenfunctions we have M{fk) = 
M(9{lk) = 0'{lk) = f'k- Since the fk are bounded functions, the multiplicative property of 
M gives M{fk, ■ . . fk„) = f'k, ■ ■ ■ f'k^ for all /ci, . . . , A;^ G 5. Then the definition of the phase 



form (|20[) and (21) shows that both processes have the same phase form. □ 



Corollary 9.2. Let u he a positive symmetric pure point measure on G and Z the associated 
group of relators. If {Zn) = Z, then two stationary spatial processes with diffraction oj are 
spatially isomorphic if and only if the first n moments of their phase forms agree. 
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Proof. This is a direct consequence of the previous theorem and Lemma 7.10 □ 



Remark 9.3. (a) The above theorem deals with isomorphism between two processes. We can 
also characterize the automorphisms of a given pure point stationary spatial point processes 
with Bragg peaks S and eigenfunctions fk, k S S. These automorphisms are in one to one 
correspondence with elements from T in the following way: Any character cj) £ T yields an 
automorphism of {N, X, /x, T) via fk i— ?• 4'{k)fk- Conversely any automorphism of {N, X, fi, T) 
must necessarily map map fk to a multiple 4'{k)fk with (j){k) £ U{1). The arising function (j) 



must then belong to T as (20) gives 4>{ki) . . . 4>{km) = 1 whenever ki + ■ ■ ■ + km = 0- 

(b) In the previous theorem, the isomorphism M between isomorphic pure point stationary 
spatial processes is not unique. However, if Mi,M2 are two such isomorphisms then M := 
M^^Mi is an automorphism of the first, say (A^, X, fi, T). By (a) of this remark, we conclude 
that the isomorphisms are unique up to elements of T. 

10. From phase forms to stationary processes with pure point spectrum: The 

torus approach 

In this section we show how to construct an ergodic full stationary process with a given 
pure point measure as its diffraction measure w and a given phase form a* as its associated 
phase form. This yields a canonical model realizing a given diffraction measure. 

We assume that we are given a a locally compact Abelian group G, a pure point positive 
symmetric backward transformable measure measure u on G. These data give rise to 

S ■.= {k£G: uj{k) > 0} 

and £ = (5)group C G. There is an associated group of relators Z as discussed in Section [Tj 
We are also given a* £ Z. 

We are going to construct a stationary process corresponding to (w, a*). 



From Prop. 7.5 we know that there is an elementary phase form a : J-{S) — > U{1) (i.e. a 
character of J~{S)) which restricts to a*: in other words, we can find a so that 

a*{ki,..., km) = a{ki) ■ ■ ■ a{km) 

whenever [k] = 0. (This a is not unique, but the ratio of any two of such elementary phase 
forms is a character on E^-) Thus, we can assume without loss of generality that a* comes 

from an elementary phase form a G ^{S). Since uo is backward transformable there is a 
measure 7 on G so that for all F G Cc(G) we have 

(22) l\F\^du= / F*Fd7<oo. 

Jg Jg 

We consider £ to have the induced topology from G and, as usual, let £d denote the same 
group £ but with the discrete topology. The dual of £d is a compact Abelian group T. We 
denote the Haar measure of total volume 1 on T by If and the counting measure on £ii by Id- 



Then, (22) can be reformulated as saying that 

{k^ F{k)u{ky/^) e£\£). 
For F € Cc(G), define a function na{F) on Ed by 

na{F) = Y,F{k)a{k)u{k)'/Hk 
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where we take the positive square roots and Ifc is the function on £fi whose value at /c is 1 
and which takes the value everywhere else. We have 



fce5 k&S 

= l\F\'^du)= / F*Fd7<oo, 
JG Jg 

which implies in particular that na{F) G L?'{£di Id)- The Fourier transform provides a funda- 
mental isomorphism between L'^(£ij_,lij_) and L'^{T,lj) taking to the character defined by 
k on T. We usually denote this character by Xk- Thus we obtain, by applying the Fourier 
transform to na{F), 

(23) NaiF) := F{k)a{k)u{kYl\k € ^^(T, h) 

kes 

with 



(24) {Na{F)\Na{F)) = / Na{F)Na{F)dlj= / n,(F)na(F) d/^ = / F * F dj . 

Jt JSd JG 

This provides us with the mapping 

and shows that for any compact subset K C G we have ||A''a(F)||2 < (|7|(i^ + i^))^/^| |-F| |sup 
for all F G C]?(G) with support inside K, which shows that Na is continuous. 

The continuous embedding of — )■ G defined by inclusion leads to a dense homomorphism 
G ^ T. Then each function on T "restricts" to one on G and in particular we have an obvious 
meaning to the functions 1 1— t- Xkit) for all i G G and for all k G S^- We also obtain a natural 
ergodic action of G on T. For t G G and ^ G T, 

t.^:k^ Xkimk) 

for all k G Sd- For F G Cc(G) and k e £d 



F{k)xk{t)= / F{x)xk{x)xk{t)dlG{x)= / {TtF){x)xk{x)dlG{x) = TtF{k). 
Jg Jg 

(Here we use the invariance of the measure Iq.) 

This leads by a straightforward calculation that TtNa{F) = Na{TtF), which shows that Na 
is a G-equivariant map. 

Thus given (uj,a*) we obtain an ergodic spatial stationary process {Na,T,lj). The above 
considerations also show that uj is the difi^raction measure associated to Na- 

We will show next that Na is full and has associated phase form a*. To do so we consider 
the map 

e : L^{G,u) L^{T,h),e{h) ^Yh{k)a{k)uj{kf/'^Xk- 

kes 

Then, 9 is an isometry with 9{F) = Na{F) for all F G Cc(G). By uniqueness of the dynamics- 
to-diffraction map we infer that 9 is the dynamics-to-diffraction map. This gives in particular, 
that 9{lk) = a{k)uj(kY^'^Xk for all k ^ S. These are the functions fk of the previous sections. 
From this equality and the fact that a is a character on J-{S) it follows easily that a* is the 
phase form associated to Na (see (20) and (21)). Moreover, Proposition 3.8 now shows that 
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a{k)uj{k)^^'^Xk = ^(Ifc) belongs to the closure of the linear span of the N{F), F G Cc(G), for 
any k £ S. As products of these Xk form an orthonormal basis of -L^(T, Ij) we obtain fullness 
of the stationary process Na- We summarize the conclusions of this section. 

Proposition 10.1. Let a pure point positive symmetric backward transformable measure oo on 
G be given and a* be a phase form. Then for any choice of elementary phase form a restricting 
to a*, (Na,T,lj,T) is an ergodic full stationary process on G with diffraction measure to and 
phase form a* . 



The pure point process constructed in Prop. 10.1 depends on the choice of the elementary 
phase form a that we choose to represent the phase form a* . What happens if w e ch oose 
another representing elementary phase form, b G J-{S)1 We know already by Thm. 



9.1 



that 

the resulting process will be isomorphic, and since both use the same constructed dynamical 
system X = T, this difference between the two processes is in effect an automorphism. Also, 



Proposition 7.5 (ii) says that the ratio u = b/a is an element of T. From the perspective of 
Na and N^^ hten is in effect a translation of Na ba u S T i.e. 

(25) Nb{F) = Y.F{k)b{k)u:^'\k = Y.F{k)a{k)u:'l\k{u)xk 

kes keS 



Y,F{k)a{k)u^/^T.u{Xk) = {T^uNa){F) . 
kes 



Remark 10.2. If u = b/a above is actually continuous with respect to the original topology 
on £ then n, being a character on £, lifts to a character on £ and then by general character 



theory to all of G. Thus u can be identified with an element of G. Now (25) becomes 

(26) Nb{F) = ^F(fc)6(fc)a;i/2^fc = ^F(A;)a(fc)a;V2^fc(n)xfc 

fce5 keS 

= ^ f^{k)a{k)uj^/\k = Na{T^uF) . 

k£S 

This time the difference between Na and Nb is a translation on G. This situation occurs in 
the periodic situation, when G is compact and G is discrete. 

Corollary 10.3. Let uj be a positive backward transformable pure point measure on G and 
Z the associated group of relators. If {Zn) = Z, then two mappings a,b : S — > U{1) define 
isomorphic stationary pure point processes if and only if their first through nth moments are 
equal. 



Proof. This is a direct consequence of Corollary |9.2[ □ 

11. The homometry problem 

In this section we discuss a main result of the paper viz our solution to the homometry 
problem for pure point diffraction. 

Definition 11.1. Given G and a measure u on G, we let AA(G, w) denote the set of all spatial 
isomorphism classes of ergodic full stationary spatial processes with diffraction u satisfying 



Assumption 6.4 
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The considerations of the previous sections rather directly prove the following result, which 
is in effect a solution of the homometry problem for pure point diffraction. 

Theorem 11.2. Let G be a locally compact abelian group. Let u be a positive backward 
transformable pure point measure on G and Z = 2(u;) the associated group of relators. Then, 
the map 

Z ^M{G,uj), 
a* ^ [{Na,T,h,T)], 
where a is any elementary phase form representing a* , is a bijection. 



Proof. By Proposition 10.1[ {Na, T, Ij, T) is an ergodic spatial stationary process with diffrac- 



tion u and phase a* . Moreover, by Theorem 9.1, the map is well-defined (i.e. independent of 



the choice of the elementary phase form representing a*) and one-to-one. It is surjective by 



Proposition |8. 1| and Theorem 9.1 □ 



We note the following immediate consequence of the previous theorem and Corollary 9.2 



(see Corollary 10.3 as well). 



Corollary 11.3. Assume the situation of the theorem. If (Zn) = Z, then two ergodic spatial 
stationary processes over G with diffraction measure to are spatially isomorphic if and only if 
their phase forms have the same m-th moments for m = 1, . . . ,n. 

Remark 11.4. Although it is obvious from the theory of pure point spatial processes that 
we have developed, it is worthwhile noting explicitly that the problem of classification of pure 
point spatial processes with a given uj depends on only the positions of the Bragg peaks, 
namely S. Once one has S, one has the group Z, which depends only on the relationship of 
S to the subgroup of G that it generates. And in fact this purely algebraic problem is in 



reality the entire homometry problem! See in § 12 2 for a ID case when G is finite. 



12. Compact Groups 

Cases when G is compact (or even finite) are important since they arise in the study of 
periodic and limit periodic cases of diffractive structures, and their relative simplicity allows 
us to see more clearly the nature of the spatial processes that we have introduced. In this 
section we first look at what happens in general when we are given that the group G is 
compact. 

We then move to some examples. These examples revolve around the case that G is finite 
and/or around one of the most basic of all diffraction measures, namely the diffraction of 
Z, 6z. We also look at an interesting result of Griinbaum and Moore that involves rational 
periodic diffraction on the line. The fact that there is anything to say in this remarkably 
simple case and that it seems to be very difficult to extend their results to a two dimensional 
setting is a bit of a testimony to the difficulties that are inherent in the homometry problem. 



12.1. The compact setting. We now look at the method of §10| in this compact case. We 
assume that uj is positive, pure point, centrally symmetric (a;(— A;) = uj{k), for all k £ S), 
and backward transformable. Let S be the support of uj. It is convenient to assume that S 
generates G as a group, for otherwise we may use (5) and (5) instead. Then £ = G = £d and 
the carrying space for our spatial process is 

X = T ■.= £d = G. 



38 DANIEL LENZ AND ROBERT V. MOODY 

When we wish to distinguish two roles of G, first as a group of translations and second as 
the set of states of dynamical system, we shall use notation like t and ^ = respectively. 
Moreover, for any character /c on T (i.e. k G G) we will denote the map 

by k as well. 

When we come to the transformability into the autocorrelation measure 7, we require that 
for all F £ CciG) = C(G), 

(27) 7(F *F)= f = V oj{k)\F\'^{k) = V oj{k)F7p{k). 

Jg ^ 
fcGG fceG 

Thus the autocorrelation Fourier dual to lij is given by 

(28) -i ■=^uj{k) k = ^u:{k) k 

provided that this is indeed a measure. Since G is compact, and 7 is supposed to be a Borel 
measure, this measure must be finite. 

Assuming transformability we can, for any character a : J-{S^ — > f^(l)) form the associated 
stationary process 

Na : Cc(G) ^l2(X,//), 

namely, 

NaiF) = Y^F{k)a{k)uj^/^{k)k 

for all F G C{G). Note that Na{F) belongs indeed to L'^{X,fi) due to ([27]). 

Our interpretation is that Na represents some sort of density on the space G. Each point 
G X represents an instance of the yet unspecified structure which can be paired with 

elements of Cc{G) to give 

(29) {^t,F) = Na{F){^t) = Na{F){t), 

where the left hand side needs to be given an interpretation. It turns out that such an 
interpretation can be given easily if we make the finiteness assumption that 

(30) a;GL^(G,/g). 

Note that this assumption immediately implies that uj is backward transformable as it yields 



that 7 given in (28) is indeed a measure. 

Given (30), we can define the function pa G L'^{X) by 

(31) pa:=Y,a{kMky/'k 

fcGG 

by standard theory of Fourier series. Then, a short calculation invoking unitarity of the 



Fourier transform (see e.g. (23) as well) gives 

Na{F){t) = Y.J{k)a{kMkf'\Kt) 
fceG 

(32) = / TtPa{x)F{x)dlG{x). 
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Note that we have indeed pointwise (in t) existence of Na{F)(t) as both F and u are square 
summable i.e. belong to £'^{£ii)- Comparing with (29), we find that we can identify with 
the function TfPa and the pairing between and Cc(G) with ordinary integration. 
When can this 'density' be interpreted as a measure pa on X = G? 

The idea is that for some £ X (which we can take to be G G) the structure is pa and 
that as we translate around the translate TtPa represents the structure at ^j. Usually one 
would have to assume by ergodicity the denseness of this orbit, but in the compact case as 
we have it here, the orbit is the entire space X. This entails that 

(33) {TtPa,F) = Na{F){t). 

where and then TtPa = ^j^^Q{k,t)a{k)uj{k)^^'^k. This gives (29) when we treat pa as a 
measure. 

As G is compact the L^-function TtPa belongs to as well and we can consider it to be 
the measure 



Ttp, 



Thus, in this situation we can realize the process as a measure process. 

If ^a;(A;)^/^ < oo, then pa will even be a continuous function. If Y^^i^)^^"^ = °o, the 
situation becomes different. Such a situation can easily arise: it suffices to find a and co such 
that Pa does not belong to L°°{X). Let for example G be compact admitting a function 
p € L^(G) \ L°°(G). As p belongs to L^, we can expand p in a Fourier series 

p = Y^a{k)uj{k)^/'^k 
k 

with a{k) £ U{1) and uj{k) > with X^w(A;) < oo. The diffraction of p is seen directly to be 
equal to oj. 



For more on interpreting Na as some sort of measure induced density on G see ^15.3 



13. Examples around the diffraction measure uj = 5z 

In general the number of solutions to the inverse problem is vast, even for the most basic 
diffraction patterns. In this section we consider the most famous example of a diffraction 
pattern, the diffraction of the integers. More properly this is the diffraction of the Dirac 
comb 5z = Ylxez which represents a point density of one at each of the integers. Its 
diffraction is also u = dz- 

Crystallographers have long relied on additional information (like the periodicity of crystals 
and their known constituents) to find solutions to the inverse problem. We shall see that 
by imposing periodicity and/or imposing arithmetic conditions on the solutions, the inverse 
problem we can do the same. One of the advantages of our approach, that applies to all locally 
compact Abelian groups, is that it also applies to compact and to finite groups. Although our 
object of attention is 6z, much of this section is devoted to cases in which G is finite and to 
showing what the theory then looks like and tells us. We also deal with the case G = U{1). 

In the case of diffraction from a periodic set or periodic crystal, it is customary to compute 
the diffraction directly from the originating measure describing the density p, by taking the 
square-absolute value of the Fourier transform of p. In the case of 5z this approach gives the 
answer directly as a consequence of the Poisson summation formula. This also works directly 
for finite groups. 
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We start with the pair of dual groups j^'L/'L and 'L/M'L where M is a positive integer. 
The duahty is most conveniently written down in the form 



for X G jyZ/Z and k G 'L/M'L. In the context of this paper, one of these groups plays the 
role of G and the other of G. Most of the time it will be G = -jjLjL and G = L/ML. To 

help clarify things, elements of G will be written in the form Xx or Xk accordingly. 

The idea here is to use the fact that the diffraction 5i is periodic and then represent it in 
L/ML and look for all solutions for an originating distribution p in jjL/L. The resulting 
distribution will then be interpreted as a periodic (modulo L) measure on j^L. In this way we 
pick up periodic solutions to the inverse problem. Alternatively we could do things the other 
way around and consider 5^ represented in j^L/L and find solutions to the inverse problem 
in L/ML which we could interpret as periodic solutions in L. However, we shall see that does 
not lead to new solutions to the inverse problem. 

We note that for finite G and G the Haar measure we want to use is counting measure 
normalized to total measure 1 by dividing by the order of the group. For Xk S L/ML, simple 
calculation gives the expected result that Xk = ^k (note that Xk is a function on G and its 
Fourier transform is a function on G). 

Since the groups are finite, there is no complication in computing diff'raction. For a density 
or distribution of density p on G, the autocorrelation is 7 = /3 * p and the diffraction is 

u: = {p*p)^ = 

Remark 13.1. We should note that in general the approach via the square- absolute value 
of the Fourier transform is not applicable to aperiodic structures. As A. Hof has shown [19j 
and [20] , in the study of aperiodic crystals one usually meets distributions of density that are 
not Fourier transformable. 

Our approach is to follow the lines laid out above and for given diffraction a; and elementary 
phase form a define the associated density via (|31[): 



This is a finite sum and there are no questions about convergence or ambiguities about how 
to interpret it. 

13.1. Specific values of M. 

13.1.1. The case M = 1. This seems too trivial to consider, but it is interesting nonetheless. 
G = \L/L, G = L/L, u) = (5o. Here 5 = {0} G L/L and T is generated by e(0) = 1. So J" 
is the trivial group and its only character is the trivial character. Also Z = T . The density 
measure p defined by the trivial character a is 0(0)07(0)^/^x0 = Xo which is simply the identity 
function on G. Thus we obtain the expected solution to the inverse problem: Xo = i^o whose 
absolute square is again uj = 60. The function xo is the identity function on jL/L, which is 
also the measure on G. 



{k, x) = e 



(34) 




fceG 
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13.1.2. The case M = 2. (1) G = ^Z/Z, G = Z/2Z, oj = Sq + Si. The two elements of 

G are xo stnd xi- We have S = {0. 1} = G and J- is generated bv e(0) = 1 and e(l), with 
e(l)2 = 1 since —1 = 1 mod 2Z. Thus J-" is cychc of order 2. There are only two characters 
on J^, i.e. elementary phase forms, 0+ : e(l) 1— t- 1 and a_ : e(l) —1. The exact sequence 

l^Z^T^£ = G^l, 

comes from the mapping e(k) k which here reads e(0) H> and e(l) 1— t- 1. Thus the kernel 
Z is trivial. The significance of this is that there is only one phase form (trivial!) and thus 
the two solutions we derive from the two elementary phase forms will be the same up to 
translation. In detail, we have, using a+ 

p+ = a+(0)a;(0)V2^ + a+(l)c^(l)V2^ = Xo + Xi ■ 

Note that this is a function on G = |Z/Z. Its values are p+{0) = 2,p^{l/2) = 0. This is a 
nice transformable measure and its Fourier transform is the measure is 5o + Si on G = Z/2Z. 
Its absolute square is w, as we expected. 
Similarly, using a_ we obtain 

p_ = a^{0)uj{0)^/^xo + a_ (1)^^(1)1/2x1 = Xo - Xi • 

Its values are P-(0) = 0,p_(l/2) = 2, so the density distribution p_ on G is just translation 
by 1/2 of p-i- found in the case of a+. The Fourier transform is Sq — Si and again the absolute 
square of this is oj. 

(2) Next we reverse the roles of G and G: G = Z/2Z, G = ^Z/Z, lo = Sq. Notice that 

G represents the reduction of the integers and the half-integers modulo Z. The reduction of 

the diff'raction 5z modulo Z is thus supported only on the Z part which is represented by the 
class of modulo Z. Following the general discussion of compact groups we should, at this 
point, replace G by the subgroup (5) = {0}. This then reduces us to the case M = 1 already 
discussed. 

(3) It is interesting to consider G = Z/2Z and G = |Z/Z with the diffraction 00 = 
Sq + S1/2. This may be considered as the reduction modulo Z of (5z + (5i .j.^. The two elements 

of G are xo and Xi/2- The discussion parallels that in (1). We have S = {0,1/2} and T 
is generated by e(0) = 1 and e(l/2) with e(l/2)2 = 1 since -1/2 = 1/2 mod Z. Thus 
J- is cyclic of order 2. There are only two characters on J^, i.e. elementary phase forms, 
a+0 : e(l/2) i-)- 1 and a_ : e(l/2) 1— )• —1. The exact sequence 

1 — > Z — >T — >£ = G — ^1, 

comes from the mappinng e{k) 1-^ k which here reads e(0) ^ and e(l/2) ^ 1/2. Thus again 
the kernel Z is trivial and the two solutions we derive will be the same up to translation. 
Using a+ 

p+ = a+(0)a;(0)V2^ + a+{l/2)uj{l/2f/''xyi = xo + X1/2 • 

Thus the distribution on Z/2Z is p+(0) = 2,/9+(l) = 0. Its Fourier transform is the measure 
5q + (5i/2) whose absolute square is oj. 

Going back to Z we have 2821 on Z whose diffraction is |(2(5ig) = 8z + Si_^^, which agrees 
with what we obtained after removing all periodicity. 
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Similarly, using a_ we obtain 

p_ = a_(0)a;(0)V2^ + a.{l/2)u{l/2Y'^x^2 = Xo - X1/2 , 

and the distribution on Z/2Z is p-(0) = 0, p-(l) = 2. Then 'pZ = 6q — ()i/2- Taking absolute 
squares we again obtain the diffraction u. 

Again going back to Z, this time we have 2821+1 on Z whose diffraction is 5z — <^i+2' which 
agrees with what we obtained after removing all periodicity. 

13.1.3. The case M = 3. Here something more interesting happens. G = ^Z/Z, G = 
Z/3Z,a; = 5o + '^i +<^2- The three elements of G are xo; Xi) X2) with x/c(i/3) = {k,j/3) = Cs'^) 
where C3 = e^''*/^. We have S = {0, 1, 2} = G and T is generated by e(0) = 1, e(l), e(2), with 
e(2) = e(l)~"^ since —1 = 2 mod 3Z. Thus ~ Z with e(l) as a (multiplicative) generator. 
There are infinitely many characters on J^, i.e. elementary phase forms, one for each u & U(l): 
Gu ■ e(l)" G C/(l) and e(2)" h+ n"" for all n G Z. The exact sequence 

1 — >Z — > T £ = G ^ 1, 

comes from the homomorphism that maps e(l)'^ k mod 3 (which includes 6(1)"*^ = 
e{-l f ^ -k mod 3). The kernel Z ~ 3Z via the isomorphism ^ ~ Z. 
For each u e U{1) we obtain 

(35) pu = a,(0)w(0)V2^ + a„(l)w(l)V2^ + a„(2)a;(2)V2^ 

= Ixo + +^^Xi • 

Thus the distribution on ^Z/Z is Pu{0) = 1 + u + u, pu(l/3) = 1 + n("| + n^3, and pu(2/3) = 
1 + ^^Cs + ^Cl • Since Xfc = ^fc we obtain 'p^ = 60 + U62 + Its absolute square is w as 
required. 

We note that replacing u by or changes the distribution p by translation, which 
corresponds to the fact that the elementary phase forms map 3 : 1 to phase forms. 

Thus there is a one parameter family of solutions to the inverse problem for lu, which are 
supported on Z and periodic of period 3. 

13.1.4. Other cases of M . If we look at the inverse problem to + • • • + 5m-i in Z/Z 
the situation evolves in a natural way. At M = 4 the value 2 leads to e(2)2 = 1, and 
e(l) = ''(3)""^ !->■ -u G U{1) is arbitrary, so we end up with solutions Xo + ^iXi/4 =t X1/2 +^X3/4■ 
At Af = 5 we are get two free parameters u,v G U{1)^ one relating to the pair 1, 4 and one to 
the pair 2, 3 (integers taken modulo 5), and the solutions are Xo+^Xi/5+''^X2/5+^X3/5+^X4/5- 
The results for higher M follow a similar pattern. 

13.2. What happens when G = U{1),G = 7/1 We again consider the diffraction but 
now use Z as G and C/(l) as G. Then iS = Z and T is generated by {e{k) : k & Z} with the 
single set of relations e(— A;) = e{k) (along with e(0) = 1). Then we have elementary phase 
forms a : T — > [/(I) with a{k) G C/(l) arbitrary except that a(0) = l,a(— A;) = a(A;). Each 
such a leads to a formal solution 

Pa = Y^ a{kMk)'/^Xk = E «(A;)e-'"*'^-" , 
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This is simply a formal Fourier series on Z. However, given that for all k, \a{k)\ = 1, there is 
no convergence as a Fourier series. This is expected since in the simplest case when a{k) = 1 
for all k its Fourier transform is the unbounded measure Sz- 
However, the solution for A^^: 

Na{F) = F{k)a{k)oj{kfl'^ ^ = Y^ F{k)a{k) k 

for all F € C{U{1)) is well-defined as an L^-function, and we can learn something from it. In 
the case that a is the trivial elementary phase form we obtain 

Na{F)it) = YF{k)ik,t) = Fi-t) = 5oiTtF) . 



Thus the associated density is ^o, which is what we expect. Recall [12.1 here. The compact 
group T, which is U{1) here, is viewed both as a dynamical system with states S,t and as a 
group of translations t £ U{1) which translate the states about. Na{F){t) is something that 
gives the effect of the process on the function F at the state ^j. 

Now let K = —K be a finite symmetric subset of Z\{0} and let a be an elementary phase 
form which satisfies a{k) = 1 for all k ^ K and a{k) = a{—h) G U{1) arbitrary for all 
A; e K n Z+. Now for F G C{U{1)) we have 

Na[F){t) = Y,F{k){k,t) + Y,^a{k) - l)F{k){k,t) , 
feez keK 



for all t £ U{1). Following the same reasoning as in (32) we arrive at Na{F){t) = (5o(T(F) + 
piTtF) where 

P = Y{a{k) - l)xk. 

k&K 

Remarkably the density 6o + Y^keKi'^(^) ~ ^)Xk on U{1) is indeed a solution to the inverse 
problem for 6q, as can be seen by directly computing its autocorrelation and taking its Fourier 
transform. This shows the ubiquity of quite reasonable looking solutions to the problem. 

13.3. A simple periodic case with imposed arithmetic conditions. In this section we 
consider the very simple situation of the diffraction from the integers when they are weighted 
periodically with period M by a set wq, . . . , wm-i of real numbers. In this case the diffraction 
is also periodic, and we will assume that the ambient space is G = Z/MZ and the diffraction 
occurs in its dual, G ~ -p^Z/Z, say 

(36) uj = oj{k)5k ■ 

k&IM'L 

As we noted above, making the assumption that G = TLjMTL is not entirely innocent since 
it puts considerable limitations on the kinds of solutions available - namely that they are 
periodic of period M and the distribution of intensity is entirely confined to integers, neither 
of which is mandatory. Taking a larger ambient space will in general produce more solutions 
(many more!). 



Equation (31) defines the density on G arising from uj and the elementary phase form a. 
Since its diffraction is to be w, this leads to the M equations 

(37) paix) = a{k)uj{kf''^{k,x) = w^, x£ Z/MZ 
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where S C G = -^Z/Z is the set of non-vanishing points of uj and {k,x) = exp{2Trikx) = 
^Mfcx g Q[Cj\,/], where Cm = exp27ri/M (the values of k are of the form ko/M mod Z where 
ko G Z). Notice that for finite groups, the situation we are in here, the densities pa are 
functions and also measures, since the distinction between them disappears. 

Griinbaum and Moore [IB] consider the case in which pa is rational valued, i.e. all of the 



weights Wx are rational numbers. In this case, solving the system of equations (37) for the 'un- 
knowns' a{k)u>{k)^/'^ we see that they must all lie in Q[Ca/]) and then a{k)u>{k)^^'^a{k)uj{k)^/'^ = 
uj{k) £ QICm] also. 

Let H := Gal(Q[CM]/Q) ^ (Z/MZ)^, the group of units of the ring Z/MZ. We can 
view H as both a group of automorphisms permuting the M roots of unity, or as a group of 
automorphisms of -^Z/Z; namely s G (Z/MZ)^ acts as : exp(27riA;) i— )• exp(27rzsA;) and as 
Og : k sk respectively. In the rational case we are now in, it is shown in pTS] that a{S) = S 
for all a G H. We show next why this happens: We take a(k)uj{ky^'^ = if A; ^ 5. 

If /ci G G then 

(38) ^ E ki{x)pa{x) = E «(^)^(^)'^'^ E = a(fei)a;(A;i)V2 , 

which is zero if fci ^ S. 

Let a = Us ^ H . Since pa is rational valued, pa = a{pa) = S^g^ a(a(A;)a;(/c)^/^)a(A;). Ap- 



plying the process of summation of equation (38) to the element a{ki) where ki € S we obtain 
a{a{ki))uj{a{ki)). But applying a dircectly to (38) we obtain X^^gg a(fci(x))a(/?a(a^)) = 



a{a{ki)u}{ki)^/'^) . These two are equal, and this shows, in particular, that a{a{ki))uj{a{ki))^/'^ / 
0, and hence that a{ki) G S, which is what we wished to prove. In short, if we identify S as 
a subset of -p^Z/Z, then k G S and s G (Z/AfZ)^ implies sk G S. In effect 5 is a union of 
orbits of jjZ/Z under the action of the multiplicative group (Z/MZ)^. 

This condition appears in an equivalent form in the following key result of [18j. For /i G Z, 
let ord(/i) denote its order as a group element modulo M (under addition). 



Proposition 13.2. // the density pa of (37) is rational valued then for all k G S and for all 

j €z Z with gcd(j, ord(M/c)) = 1 one has jk G S. In this case any mapping a : S — > U{1) 
satisfying the m-moment condition for all m < 6 if M is even (respectively m < A if M is 
odd) extends to a character on jjZ/Z. □ 



Using Proposition 7.6, and Corollary 11.3 and the corresponding notation, we obtain 



Corollary 13.3. Under the conditions of Proposition 13.^ the set of distributions pa whose 



diffraction is uj of (36) is classified by the first 6 moments (resp. first 4 moments) if M is 



even (respectively if M is odd). □ 

We do not know of any extensions of this result into higher dimensional periodic situations, 
and indeed it seems very difficult to do. 



Example 13.4. Here we give an example based on an example given in (18] (see also [T2 
We begin with G := Z/6Z and the two weighted Dirac combs 

(39) ll(5o + 255i + 42^2 + 45^3 + 31(54 + 14^5 

10<5o + 17(5i + 35(52 + 46^3 + 39(^4 + 21^5 
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on G, which represent two density distributions for which we would hke to find ah possible 
homometric density distributions. What is remarkable is that not only are these two density 
distributions homometric, they actually have all of their first 5 moments in common. 

Let us begin with the first of these two densities. The Fourier transform of this density is a 
function on G ~ gZ/Z. The latter is most conveniently thought of as the 6th roots of unity. 
We write the Fourier transform in the form 

(40) i5 + 25w + 42w'^ + 4:5w^ + 31w'^ + Uw^ + 6w'^) 

= {w + l){w'^ +W + l){2w'^ + 5){3w + 1) =: P{w) 

where w varies over the 6th roots of 1, Uq := {e'^'^^^/^ ; j = 0, ...,5}. Note that = 
1, so the coefficient of is 11/6 as it is supposed to be. The diffraction at A; G G is 
uj{k) = P{w^)P{w^), which works out to give the following values for oj{ky/^ : 28, y^247/3, 
,0 ,0 , a/247/3 at A: = 0, 1/6, . . . , 5/6. Thus in this example S = {0,1/6,5/6}. In fact 
Griinbaum and Moore constructed the polynomial P precisely to make the Bragg spectrum 
have extinctions at 2/6, 3/6,4/6. In the rest of this subsection we shall make the notation a 
little simpler by suppressing the ubiquitous denominators 6 which occur in the values of k. 
Alternatively we multiply everything by 6 and work in Z/6Z instead of gZ/Z. 

In any m-tuple k = {ki, . . . , km) in 5™, we can, using the equivalence relation ~, drop any 
occurrences of and drop any pair 1 , 5 that occurs amongst the kj . This leads us to conclude 
the only equivalence classes of S are of the form 0~, (1, . . . , 1)~, (5, . . . , 5)~, where there are 
any positive integer m entries in the two latter types. Since (1, . . . , 1)~ and (5, . . . , 5)~ with 
m entries are inverses of each other in S, we conclude that S~ ~ Z. Since the mapping 



4) : S~ — > £d — Ue, we have if : — > Uq (Prop. 7.3). 



Using the identification of S~ and J^{S) we have J^(5) ~ Z and we can put together the 



exact sequence (17): 



1 ^ 6Z ~ Z ^ J^{S) ~ Z ^ C/e ~ = G ^ 1. 
Dualization then gives the exact sequence (flsl) 



1< — f ~ U{l)/U6 ~ ^Z/Z i — IXS) ~ U{1) i — Z/6Z ~ T ^ — 1. 

Each element of U (1) can be viewed as a character of gZ, and when restricted to Z characters 
a(-) and e^'^*-^^'-'a(-), j = 0, . . . , 5/6, are the same. 
We now write down the densities 



fceS 



"s" a(k)u:{kfl'^k{x) = 28 + 7247/3 a(l)e-2-W6 + ^24773 a(5)e2"^^/*' 



Here x G G ~ Z/6Z. The character a of -F(5) is a character of Z and hence is an element 
of [/ (1). Recall that T{S) is generated by the tuples (1,1,...) and their inverses (5,5,..., 5) 
where 1,5 represent classes modulo 6 in Z. There are no further relations. Thus a(A;) here 
really means a((A;)~), that is, the value of a on the the equivalence class of (k) under ~. Thus 
by a(l) we mean a((l)~) and by a(5) we mean a((5)~) = a(l). There is no restriction on 
what a(l) might be in C/(l), so simply writing it as a we may write p as 

Pa{x) = 28 + V247/3ae-2™/6 + ^247/3 cZe^^^^/^ 
where a G f7(l) is arbitrary. Each value of a gives a solution to the diffraction problem. 
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Figure 1. The graphs of Pa{x) as a varies for each of the six values 
0, 1, 2, 3, 4, 5 mod 6 of x. Here a varies over the unit circle {e^'^** : t G [0, 1)}, 
indicated here by the interval [0, 1) and the six graphs give the corresponding 
values of the coefficients or weights of the deltas that make up the distribution 
on Z/6Z. The original weighting distribution 11, 25, 42, 45, 31, 14 occurs when 
a ~ 0.443099. The second density of ([39]) with coefficients 10, 17, 35, 46, 39, 21 
also occurs, when a ~ 0.520310. 



The following figure illustrates the graphs of pa{x). The situation that we began with 
occurs when a ~ 0.443099. Another solution indicated in [18.j is 10,17,35,46,39,21 which 
occurs at a ~ 0.520310. 

Now we are going to show that in this example, {Zq) = but {Z^) C Z. The argument 
is quite simple. First note that (1, 1, 1, 1, 1, 1) E Zq but cannot be reduced, so its reduced 
length is 6 (note that denominators are still being suppressed). So {Z^) / Z. Now consider 
any k = [ki, . . . , kn) G Zn with n > 6, which is of minimal length and can not be split into a 
product of two (non trivial) elements from Z. There can be no occurrences of since those 
can be discarded. Since {1, 5} and {2, 4} are annihilating pairs modulo 6, we cannot have both 
1 and 5 in k and neither can we have both 2 and 4. Neither can we have 2 appearing three 
or more times since (2,2,2) G and similarly for 4. Likewise 3 can appear at most once. 
Suppose 2 or 3 appears in k. The possibilities are 2, 22, 3, 23, 223. The remaining symbols are 
just Is or 5s. But 21111, 2211, 3111, 255, 3555 then show that k is not reduced. The argument 
with 4 is much the same. So we can assume that none of 0,2,3,4 appear and not both 1 and 
5 appear. This leaves just (1, . . . , 1) and (5, . . . , 5) as the possibilities for reduced words, and 
in either case the number of entries must be a multiple of 6. 

The consequence of this is that any phase form is uniquely determined by its first through 
sixth moments. For if two have these moments in common then taking corresponding elemen- 
tary phase forms a,b we find that ba~^ £ J'i'S) and satisfies the first 6 moment conditions. 
Since (Zq) = Z, ba~^ lies in T and so the corresponding phase forms were equal. 

It is a surprising fact, pointed out in [18], that the two densities of (39), IISq + 255i + 
42(52 + 45(53 + 31(54 + U65 and W60 + 17(5i + 35(52 + 46(53 + 39(54 + 2I65 indicated above, have 
all of their first through fifth moments in common, and it is only at the sixth moments that 
they look different. 
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14. An alternative construction of spatial processes from pure point 

diffraction 



In ^ 10 we offered a direct method of constructing a spatial process from the original data 
consisting of a positive backward transformable point measure w on G and a phase form 
a* (or equivalently an elementary phase form a). This approach can in a certain sense be 
understood as an elaboration of the theorem of Halmos - von Neumann which tells us in 
advance that since the required process will have an associated dynamical system with pure 
point spectrum, it can be modelled (measure theoretically) on a compact Abelian group. 

In f|5] we have used the Gelfand method to produce pure point spatial processes in the 
splitting of spatial processes into their pure and continuous parts. 

In this present section we sketch out how we can also use the Gelfand method to produce 
a pure point spatial process out of the given pure point measure w and the elementary phase 
form a. The method ultimately reverses many of the steps that we saw earlier in which, 
starting from a process we constructed the diffraction and the diffraction-to-dynamics map. 
For this reason we do not give proofs but rather provide an outline of intermediate steps. 



We use the notation of ^ For simplicity of presentation we shall assume that G 5. 
Define c : S — > C by 

c{ki, ...,kn)= a{ki, . . .,kn)uj{ki)^ . ..uj{kn)^ , 

or more compactly, c(k) = a(k)a;(k)2, where for each square root we take the non-negative 
root. We define c(0) = 1. Note that c(k(-k)) = a;(k). 

Let C[Fk : G 5] be the free associative algebra in the variables F^. For elements k = 
{ki, . . . , km) where the kj G S, we will write Fk := Fk^ . . . Fk^ . We define an algebra <f -grading 
on C[Ffc : /c G 5] by assigning degree k to each F^, so Fk gets degree [k]. 

Let / = I{a) denote the ideal of C[Ffc : /c € 5] generated by all the elements of the 
form Fk — c(k) for all k = {ki, . . . km) S Z, i.e. for all k with [k] = (this includes Fq — 
ioiO)^/"^!, FkF^k - ^{k)) and let 

V = V{a) := C[Fk : k G S]/I . 

We let 

a : C[Fk :keS]^V 
be the natural homomorphism and set 

fk := a(Ffc), /k := a(Fk) = fki ■ ■ ■ fk^ 

for k = /ci + • • • + km- 

Since c(/ci, . . . , km) is independent of the order of the kj, we see that V is commutative. 
Furthermore, since the relations imposed on C[fk ■ k £ S] are homogeneous of degree 0, the 
resulting algebra V{a) is also graded by £. We denote its space of elements of degree k £ £ 

The space C[/fe : /c G 5]^ of elements of degree k G is the linear span of all the elements 
/k with [k] = K. If [k] = [1] = K then k(— 1) G Z and the relations defining I lead to 

(41) /k/_i = A, . . . fU-i. ■ ■ . f-i„ = a(k(-l))a;(k)V2^(-l)i/2 g 1 . 

From this we see first that = C/^ where is any one of the elements /k with [k] = k. 
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We shall always assume that we are working with a phase form a for which V{a) is not 
reduced to the trivial ring {0}. Such cycles always exist, for example the trivial cycle, which 
is 1 everywhere. 

Since the algebra is f-graded, non-zero elements of different degrees are linearly indepen- 
dent and 

p = p(a) = 0C/.. 

Define an conjugate linear form on C[fk ■ k £ S] by sesqui-linear extension (with conjugation 
on the second variable) of 

" c((k)(-l))if [k] = [l]. 



(/k I /l) — {fkl ■ ■ ■ fkm \ fh ■ ■ ■ fin) 



otherwise . 



It is a straightforward exercise to show the identity 



(/ U) = (5 I /) 
for all f,g€ C[fk : k £ S]. We shall write ||/|| for (/ | f^^^. 

Proposition 14.1. Let h,l,m be tuples of elements ofS. 

(i) (/k|/k) = c(k(-k))=a;(k)>0; 

(ii) (1 I 1) = 1; (/i I 1) = c(l) when [1] = 0; 

(iii) (/k/m I /l) = (/k I /-m/l); 

(iv) (/k - c(k)l I /i) = when [k] = 0; 

(v) The kernel of {■ \ ■) is the ideal of relations that we have imposed on C[fk ■ k £ S]; 

(vi) (• I •) induces an inner product (for which we use the same notation) on V. 

Now V is an inner product space. We wish next to consider V acting on itself by multipli- 
cation and to make it into a normed algebra by use of the operator norm. The details are as 
follows. A typical element of V can be written as a finite sum x = b\f\ with h\ S C, where 
there is only one summand of any particular degree (i.e. [1] ^ [1'] if 1 7^ 1'). Then we have 

ini^ = En'ii/i||' = Envi). 

Let k = (/ci, . . . , km) G 5"". Then /kX = Yj ^i/k/i = ^i/ki with degrees [k] + [1] and 
ll/kxiP = (5]6i/kilE^i/ki)'/' 

Thus 

(42) ||/kx|| = ||/k||||x|| 

which shows that for each tuple k, left multiplication by /k is a linear operator of norm ||/k||- 
Now we see that the entire left representation of V on itself is a representation by bounded 
linear operators, and this provides us with an operator norm v onV with i^(/k) = ||/k|| and 



— Yli \b\\^{f\) in general. In view of (42) we always have 

y[f^) = a;({k})^/2 and ||x|| < v{x) for all x£V . 
Being an operator norm, v automatically satisfies 

v{xy) < v{x) v{y) . 
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We let A denote the completion of V with respect to z^. This is a commutative Banach 
algebra. The norm will still be called v. We identify V with its image in A. Since i^{f) > \\f\\ 
for f € V, the embedding of 7^ in ^ is faithful. 

We define an G action on 7^ (a) by 

Ttfki ■ ■ ■ fk^ = {h -\ h km){t)fk^ . . . fk^ , 

or in abbreviated form Ttfk = [k](f)/k, for all t £ G. We note that Tt acts as a unitary 
transformation relative to the inner product and we obtain a unitary representation T of G 
on V. Each Tt is an automorphism of ^ as a Banach algebra. 

At this point we can apply Gelfand theory to produce the space X = Spec(^). which is 
the set of algebra homomorphisms vr : A — > C. Then X is non-empty and compact in the 
Gelfand topology This is the weakest topology on X making continuous the evaluation 
mappings ??j : vr i— )• 7r(/) for each f £ A. This means that the open sets are generated from 
the open sets 

U{f,V) := {vr G X : 7r(/) G 1/} , 

where f £ A,V open in C. This is the way in which A may be viewed as a space of continuous 
mappings on X, namely /(vr) := '/r(/), and this is the point of view that we take from now 
on. 

One has for all vr G X, / G A, |vr(/)| < (see [23, Ch. 18, Prop. 1) from which for the 
sup norm we have ||/||oo < ^(/)- 

For all tuples k from S and for all vr G X, 

|/k(7r)| = ||/k||. 

Since Tj acts as an algebra automorphism of A, T^vr = vr o T_j is another homomorphism 
of A into C, and so is in X. So G acts on X. Since Tt just permutes the defining open sets 
of the topology of X, this action is continuous. 

Lemma 14.2. For each g £ A the mapping 

G^C, u{Ttg - g) 

is uniformly continuous. 
Proposition 14.3. The mapping 

Gx X — > X 

defined by the action of G on X is continuous. {X, G) is a topological dynamical system. 

Since G acts on X, it acts naturally on C{X) too. V already has an action of G and 
acquires another one when treated as a subalgbra of C{X); however, not surprisingly, the two 
actions of G on "P the same. 

Proposition 14.4. V is dense in the space C{X) of all continuous functions under the sup- 
norm. 



The next objective is to create an ergodic invariant measure for {X, G). We define a linear 
functional n = HaOnVhy linear extension of 



c(k), if[k]=0; 
otherwise 
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for all tuples k from £. Note that this is well-defined - the relations amongst the /k for 
equal values of [k] are consistent with the definition - and one can see immediately that /i is 
G-invariant since it picks up only the 0-eigenspace of V. 

We need to extend this to a linear functional defined on all of C{X). Let be a 

van Hove sequence for G. For each / G C{X) define a new function M(/) on X by 



M(/)(7r) := hm —L- j {Ttf){7r)dt, 

n-J-oo vol An J a 



where integration is with some prefixed Haar measure on G and vol(A„) is the value of this 
measure on An. 



Lemma 14.5. For all f EV, M{f) is a constant function. Moreover, 



c(k), if [k] = 0; 
otherwise . 



This shows that jj. is the same as the mean M on V. But, as we now see, M is defined on 
all of C{X) and is a bounded linear operator on C{X): 

Lemma 14.6. For all g G C{X), ||M(5)||oo < IblU- For all g G C{X), M{g) is a constant 
function. 

In this way, M produces a measure on X, which coincides on V with defined above. We 
also call this measure /x and note the fact that we have established, namely for all g G C(X), 
and for all tt G X, 

H{g) = lim — ^ / {Ttg){'K)dt, 

n-J-oo volA„ Ja^ 

which is a form of ergodic theorem. 

Proposition 14.7. For each k e £, the space of K-eigenfuctions in C{X) (for the action of 

T on X) is one dimensional and it is spanned by any of the functions for which [k] = k. 
In particular the 0-eigenfunctions are constant functions. The complex conjugate of /k as a 
function on X is /_k = c(k(— k))/j^^. 

Proposition 14.8. /j, is G-invariant probability measure on X and {X, G, fi) is a pure point 
ergodic dynamical system. 

We now wish to define the ergodic spatial process A'^. To do this we first define the 
diffraction-to-dynamics map. For each k E S, let Ik ■ G — > C be the function which is 1 at 
k and everywhere else. Since the measure u is only non-zero on the points of S, it follows 
that every function in L'^{G,lo) can be represented in the form 

^^Xfclfe, where \xk\^oj{k) < oo . 

A;e5 keS 

The diffraction-to-dynamics embedding is the mapping 

e = ea : L^{G,uj) L^{X,n) 

^Xfclfe ^Xkfk- 

kes kes 
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Comparing these two L^-spaces we have 

i^Xklk I ^Xklk) = I XkXi Ifc liduj = ^ \xk\'^uj{k) = (Y^Xkfk I ^Xkfk) , 

'^'^ k,i k 

which shows that is an isometric embedding. 

Define an action U of G on L'^{G,uj) by setting, for each g G L^(G,cj), Ut{g) to be the 
function 

Ut{g){k) = k{t)g{k) forall/fceG. 

Under this action, Ifc G L^(G,6t;) is a fc-eigenfunction for the action and 9 becomes G- 
equi variant. 

Since a; is a positive, translation bounded, and backward transformable measure on G 
there is a unique positive definite measure 7 on G for which w = 7. This measure is also 
translation bounded (and hence automatically transformable (with Fourier transform equal 
to cj), [H]). The basic relationship between 7 and uj is expressed by ([6]) for all F G Cc(G), 
see [8], Ch. 1.4, which in particular gives F G L^(G,a;). Direct calculation shows that for 
FeCc{G), UtF = JUF. 

We define the process by 

N = Na: a(G) L\X, f,), N{F) = 9{F) . 

Also N is real, linear, continuous, real, and stationary (G-equivariant). It is also full. 

Proposition 14.9. J\f = {N, X, fi, T) is a pure point ergodic spatial process with diffraction 
UJ and associated phase form a* . 

Thus, based on the positive measure a; on G and a phase form a* we have constructed a 
pure point spacial process on G. 



15. Various concluding remarks 
In this section we collect together various small items that are of interest. 

15.1. The case that maps into A spatial stationary process N = (N, X, ^, T) is said 
to have an mth moment if N{Fi) . . . N{Fm) belongs to L^{X, ^) for any F\, . . . , F^ G Cc(G). 
In this case, the mth moment is defined to be the unique linear map jJ-"^^ on Cc(G)®- • •(8)Cc(G) 
(m factors) with 

/x('")(Fi ®...®F^):= j N{Fi) ■ ■ ■ N{F^)dii. 

It it not hard to see that the measure /i determined by its moments m G N, if all these 
moments exist and the process if full. In fact, this is more or less the definition of fullness. 
In order to have a meaningful diffraction theory our running assumption is that the second 
moment exists and is a measure. Then, also the first moment is a measure (as can easily be 
seen). 

We say that N is hounded if N{F) G L°^(X, /x) for every F G Cc(G). In this case, for every 
m G N the mth moment exists. 

We are now going to discuss how for bounded processes this notion of moment is related 
to the notion of moment of a phase form introduced above. 



/ /fcl • • • fkmdfJ- 

Jx 
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Thus, let {N, X, /i, T) be an bounded, ergodic full stationary process with pure point spec- 
trum and elementary phase form a. Then, the mth moment of and the m-moment of a 
uniquely determine each other. More precisely, the following holds in this situation. 

Lemma 15.1. (a) Let n G N and Fi, . . . ,Fn G Cc(G) be given. Then, 

kieS k„&S ''^ 

where the sums exist if taken one after the other, and for any ki, . . . , kn £ S , 

a{ki, km)uj{kif/'^ . . . w(A;^)V2 ^j^^ _^ . . . _^ /.^ = Q; 
^0 ifki + --- + kmi-Q. 

(b) For each j = l,...,n, let {Fj^^} be a sequence in Cc(G) whose Fourier transforms 
converge to Ikj in L'^{G,uj). Then, 

fk,...fk„= lim hm ... hm NiFi""'^) . . . N{F.^'^"^), 

rrti— >cxD m2— >oo rrin—^oo 

where the limits are taken in and one after the other. In particular, 

I fk,-..fkjk^= lim lim ... lim f NiFt"'^) . . . N{Fir-^)dfi . 

mi— >-oo m2— >C)0 m„— >oo 

Proof. This is essentially proven in Lemma 2.1 and Lemma 2.3 of [29j for special uniformly 
discrete point processes. The proof given there carries over to our situation. □ 

Corollary 15.2. Let two bounded spatial processes with pure point diffraction and the same 
diffraction measure u be given. Then, these processes have the same m-moment if and only if 
their associated phase forms have the same m moment. In particular, if furthermore {Zn) = 
Z, then the two processes are isomorphic if and only if their m moments agree for m = 
l,...,n. 

If every element of £ can be expressed as a sum of m or less elements of S, then (^2m+i) = 
Z, and so 2m + 1 moments suffice to determine the isomorphism class of a process, see 
Remark 17.71 



15.2. A topological situation. In this section we look at a pure point stationary process 
N : Cc(G) — > L^(X, /i) and assume that X is a topological space and all eigenfunctions 
are continuous. Here, continuity means that for each eigenvektor k, we have an continuous 
nontrivial function gi^ with 

gk{Ttx) = {k,t)gk{x) 

for sdl X £ X and t £ G. Without loss of generality we can then normalize the gk to satisfy 

gk{xo) = 1 

for some fixed xq £ X. These normalized gk will then obey 
(43) gkigk2 = aki+k2 and g^ki = Wi- 

Of course, as discussed in Section [8j the pure point process comes with a diffraction measure 
oj a set of eigenvalues £ and a phase form a. Set T := E^,. Then the considerations of Section 
[s] provide a torus system (A^'a, T, 7")- This system can be directly calculated as follows: 
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Theorem 15.3. Assume the situation outlined above. Then, 

TT : X — > T, X I— )• (fc I— )• gk{x)), 
is a surjective G-map (and hence T is a factor of X) and 

M:L\T)^L\X),g^go7r, 
is an isomorphism of point processes. 



Proof. Due to the normahzation and (43), 7r(x) belongs indeed to T. Obviously, vr is con- 
tinuous (as £ is given the discrete topology). Moreover, M can easily be seen to map the 
character k € £ oiT to the eigenfunction gi^. Thus, M is unitary (as it maps an orthonormal 
basis to an orthonormal basis). As M is an isometry, vr must have dense image. As it it a G 
map it is then onto. Moreover, a short calculation invoking (|43|) again shows that 



Migk^gki) = gkigki- 

By usual limiting arguments, this gives that M[fg) = M{f)M[g) for all f,gG L°° . 

As M maps the character k ^ £ to the eigenfunction gk, its inverse must map the 

eigenfunction gk to the character k. Combined with the usual limiting arguments this shows 
that also satisfies M~^{uv) = M~^{u)M~^{v) for bounded u,v. 

So, if we now define the A^-function At on the torus by 

Nj := M-^N 

we have indeed an isomorphism of processes. Thus, the phase form belonging to Nj must 
then be given by a as well. By uniqueness then [Na] = [Nj]. □ 

Another issue around continuity is the possible continuity of the mappings a : S — > U{1) 
in the topology on S induced from G. For these one uses this same topology to make J-'{S) 
into a topological group. We do not know if continuity here is physically relevant but we do 
note that in this case the classification results for elementary phase forms compared to actual 



phase forms change from T to G (see the discussion around and including Remark 10.2). 



15.3. Spatial processes arising from measures on G. Our underlying understanding of 
a spatial process M = {N, X, fi, T) is that it represents some sort of 'density' on the space 
G. Each point ^ G A represents an instance of the yet unspecified structure which can be 
paired with elements of Cc(G) to give 



(44) {p^,F) = N{Fm. 

Thus Cc(G) provides a set of test functions for some sort of distribution at ^ G X. Of 
course N{F) is not necessarily defined at all points of X (it is only an L^-function on A), and 
even if it is defined it is not clear whether or not could be interpreted as any particular 
kind of distribution. We do note however, that from G-invariance we have N{F){Ttp^) = 
{T^tN{F)m = N{T^tFm,so 

{Ttp^,F) = {p^,T^tF) 

for all t e G. 

An important question, and one for which we cannot give a satisfactory answer, is when 
can this 'density' be interpreted as a measure p^ on X = G? In the theory of Delone point 
sets A this actually happens; typically X is the orbit closure of A under the local topology 
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and is the process that arises from the point measures 6^ := X^^-g^ Sx for all ^ G X. Thus 

for F G C,(G), N{F){0 ■= E^eA^i^) = ih^^)- 

We have already seen one way in which we can see from general principles that measures 
arise in §12[ But this is unnecessarily restrictive since in many basic situations cj is not 



summable. For example, in § |13.2| we looked at the diffraction of the point set 6z, when it is 
modeled on U (1) as 6o. Its diffraction is lo = 6z on Z. This is not with respect to the Haar 
measure of Z (counting measure). Here X can be identified with U{1) (since it is compact) 
and for all t G U{1), F i— t- N{F){t) is just the measure 6t. So in this case there is a suitable 
measure in spite of the lack of summability. 



The exact conditions that defined by (44) should be a measure at are that N{F) 



actually be defined at ^ for all F G Cc(G) and that for each compact subset of G there is 
a constant ck so that for all F G Cc(G) with support inside K we have N{F){^) < cj^||F||oo- 
In fact this is precisely the definition of measure via the Riesz representation theorem. 
If for some ^ G X, N{F){(^) is defined, positive, and finite for each positive F £ Cc 



then the right-hand side of (44) provides a positive linear functional on Cc(G) and so defines 
a positive measure. 

15.4. Some open questions. There are numerous questions that remain to be investigated. 
Amongst these we can point out: 

• Are there simple process theoretical conditions on pure point measures lo that allow 
one to conclude that the processes atttached to w can be interpreted as measure valued 
processes? 

• Our homometry theorem shows that the stationary processes with fixed pure point 
diffraction form a group. Does this have any deeper meaning? Can we directly realize 
multiplication / inversion of stationary processes? 

• When there are moments of all orders, these moments characterize the corresponding 
processes. There are various scenarios under which not all these moments are nec- 
essary to solve the inverse problem. For instance in [lOj it is seen that real model 
sets with real internal spaces are recoverable from their second and third moments. 
In |29j it is seen that if 5 + • • • + 5 = £ [m summands) then one needs only the 
first 2m + 1 moments to solve the inverse problem. This suggests that solutions to 
the inverse problem can be organized into hierarchies depending on the number of 
moments required to specify them. 

• The sequence of moments suggests that there is some sort of cohomology theory that 
would allow one to organize the increasing information that is added as successive 
moments are included in the picture. Such a cohomology theory would be a most 
welcome addition to this study. 
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